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1. Introduction 

The concept of folding has an important place in polymer physics. Typically one 
considers the statistical model of a polymer chain made of say n identical constituents, 
and which may be folded onto itself. The entropy of such a system is obtained by count- 
ing the number of inequivalent ways of folding the chain. The combinatorial problem of 
enumerating all the compact foldings of a closed polymer chain happens to be equivalent 
to another geometrical problem, that of enumerating meanders []!]], i.e. configurations of 
a closed road crossing a river through n bridges. To our knowledge, this is still an open 
problem, which indeed has been adressed by very few authors. 

In the present paper, we study various aspects of this meander problem. In section 2, 
we define the meander problem itself, as well as a somewhat simpler semi-meander problem 
corresponding to the compact folding of an open polymer chain. We also gather in this 
section a number of data on exact enumeration of some meanders, as well as conjectural 
analytic structure for some of these data. A meander may be viewed as a particular 
gluing of two arch configurations, representing the configuration of the road respectively 
above an below the river. Section 3 is devoted to the derivation of exact results for the 
statistical distribution of arches in arch configurations. In section 4, we introduce an 
exact recursion relation for constructing all semi-meanders, for which we develop various 
approximations leading to estimates of the entropy of folding of open polymers. Section 
5 presents an alternative description of the meander problems in terms of random matrix 
models, and gathers a few consequent results. A number of sum rules and inequalities 
satisfied by meander and semi-meander numbers are displayed in section 6, together with 
exact solutions for simpler meandric problems. In section 7, we address the question of 
irreducible meanders [[J, i.e. systems of several roads crossing a river, which are interlocked 
in an irreducible way: in this case, an exact solution is derived for both irreducible meander 
and semi-meander numbers. A few technical details are gathered in appendix A. 

2. Definitions and generalities 

2.1. Meanders 

A meander is defined as follows. Let us consider an infinite straight line (river). A 
meander of order n is a closed self-avoiding connected loop (road) which intersects the line 
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(a) (b) 

Fig. 1: The M 2 = 2 meanders of order 2 (a), and the M 2 = 2 two- 
component meanders of order 2. The infinite river is represented as a hori- 
zontal line. 

through 2n points^ (bridges). A meander of order n can clearly be viewed as a compact 
folding configuration of a closed chain of 2n constituents (in one-to-one correspondence 
with the 2n bridges), by putting a hinge on each section of road between two bridges. 
Note also that the road and the river play indeed similar roles, by appropriately cutting 
the road and closing the river. 

In the following, we will study the number M n of inequivalent meanders of order n 
(by inequivalent we mean meanders which cannot be smoothly deformed into each other 

(k) 

without changing the order of the bridges) . We will also be interested in the numbers MA 

of inequivalent meanders of order n with k connected components, i.e. made of k closed 

connected non-intersecting but possibly interlocking loops, which cross the river through 

a total of 2n bridges. Note that with this last definition, M n = m£\ The M 2 = M 2 (1) = 2 

(2) 

meanders of order 2 and the M 2 = 2 two-component meanders of order 2 are depicted 
in Fig.l (a) and (b) respectively for illustration. 

2.2. Folding a strip of stamps, semi-meanders 

As mentioned in the introduction, the meander problem is equivalent to that of com- 
pact folding of a closed polymer chain. In this section, we will instead consider the case of 
an open polymer chain, moreover attached by one of its extremities. In another language, 
the problem is nothing but that of folding a strip of stamps || ||, and leads to a 
slightly different version of the meander problem, the semi-meander problem, which we 
describe now. 

The problem of folding a strip of stamps may be stated as follows. One considers 
a strip of n — 1 stamps, the first of which (say the leftmost one) is fixed, and supposedly 
attached to some support, preventing the strip from winding around the first stamp. A 
folding of the strip is a complete piling of its stamps, which preserves the (non-intersecting) 

1 The number of intersections between a loop and an infinite line is necessarily an even number. 
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Fig. 2: The 4 inequivalent foldings of a strip of 3 stamps. The fixed stamp 
is indicated by the empty circle. The other circles correspond to the edges of 
the stamps. The first stamp is fixed and attached to a support (shaded area). 

stamps and their connections, and only affects the relative positions of any two adjacent 
stamps: each stamp is folded either on top of or below the preceding one in the strip. 

The number of inequivalent ways of folding a strip of n — 1 stamps is denoted by S n . 
In Fig. 2 (a), we display the S4 = 4 inequivalent foldings of a strip of 3 stamps (stamps are 
represented in side view, by segments, and also not completely folded to clearly indicate 
the succession of folds). 




Fig. 3: The M3 = 2 and M4 = 4 semi-meanders of order 3 and 4. The 
source of the corresponding semi-infinite river is indicated by an asterisk. 

This folding problem turns out to be equivalent to a particular meander problem, 
which we will refer to as the semi— meander problem. Let us consider a half (semi- 
infinite) straight line (river) starting at a point (source). A semi-meander of order n 
is a closed self-avoiding connected loop which intersects the half-line through n pointsi 
(bridges). Let us denote by M n the number of inequivalent semi-meanders of order n. In 
analogy with the definition of meanders with k connected components, we also denote by 

~~ (k) 

Mn the number of inequivalent semi-meanders of order n with k connected components. 
As before, we have in particular Mn = M n . The M3 = 2 and M4 = 4 semi-meanders of 
order 3 and 4 are depicted in Fig. 3. 
Let us prove that 

Sn = Mn- (2-1) 

2 The number n of bridges for a semi-meander need not be even, hence n = 1, 2, 3, ... 
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(iv) 
(iii) 



Fig. 4: The mapping from semi-meanders to folded strips of stamps, (i) 
cut the leftmost arch of the semi-meander along the river, (ii) stretch the 
open circuit into a line, (iii) identify the segments of bent river with edges of 
stamps, (iv) draw the stamps according to the relative positions of crossings. 

Starting from a semi-meander of order n, let us construct a folding of the strip of n — 1 
stamps as follows. Cut the leftmost arch of the semi-meander along the river as indicated 
in Fig. 4 (i), and stretch the (now open) circuit into an infinite line (ii). The river has been 
bent in this process, but the structure of its crossings with the line is preserved (except 
that the leftmost bridge has been erased). In a third step (iii), each segment of the river 
between two crossings of the infinite line is identified with an edge of stamp. In particular, 
the last edge of the strip of stamps is associated with the segment of river between its 
source and its first bridge, whereas the first edge of the fixed stamp corresponds to the 
infinite segment of river after the last bridge. This last choice singles out the first edge of 
the first stamp in such a way that no piece of the strip can wind around it. Finally (iv), 
the stamps can be drawn between the edges, the connection being indicated by the river. 
Note that in the process one of the bridges has been erased by the initial cut, henceforth 
we are left with n — 1 bridges, hence n — 1 stamps. This gives a one to one mapping from 
semi-meanders of order n to folded strips of n — 1 stamps, thus proving (|2.1|). 

An analogous construction allows one to relate the meander number M n to that of 
foldings of a closed strip of 2n stamps. This construction is dual to the direct equivalence 
mentioned in the previous section, in the sense that the road and the river are exchanged. 
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2. 3. Numerical data and a few conjectures 

In this section, we gather numerical data for multicomponent meanders and semi- 
meanders. 

Meanders. A few preleminary remarks are in order. The number of connected compo- 
nents for a meander of order n may not exceed n, as each connected component of the 
meander must cross the river an even number of times, whereas the number of bridges is 

(k) 

2n. We display the first values of the multi-component meander numbers MA in Table I 
below (some of these numbers can be found in [I] and ||). 
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Table I: The numbers MA of inequivalent meanders of order n with k connected com- 
ponents, for 1 < k < n < 12, obtained by exact enumeration on the computer. 



The reader will check on Table I that 




(2.2) 



known as the Catalan number of order n, denoted by a 
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(n) 

Encouraged by the relative simplicity of MA , we have tried various guesses for the 
explicit form of the M^ 1 leading to the following conjecturei values of j = 0, 1, 2, 3, 4, 5 



M (n) 

n 

M (n_1) 



(2n)! 



n\(n + 1)! 

2(2n)! 
(n- 2)\(n + 2)! 
2(2n)! 



(n-3)!(n + 4)! 

M( n " 3) = — -(n 4 + 20n 3 + 107n 2 - 107n + 15) 

3(n-4)!(n + 6)P ! 



M, 



(n-4) 



2(2nV 

— Vtt^ rr(^ 6 + 39n 5 + 547n 4 + 2565n 3 - 5474n 2 + 2382n - 672) 

3(n- 5)!(n + 8)P ; 



M( n " 5) = — -(n 8 + 64n 7 + 1646n 6 

15(n-6)!(n + 10)! V 

+ 20074n 5 + 83669n 4 - 323444n 3 + 257134n 2 - 155604n + 45360) 

(2.3) 

This leads to the more general conjecture that the number ^ of (n — j)-component 

meanders of order n has the general form 



j!(n + 2j)!(n- j 



■jTj x ^- 2 H , (2.4) 



for j > 1, where P2j-2(x) is a polynomial of degree 2j — 2, with integer coefficients. The 
first few coefficients of Pij-i read 

P 2j -2(x) = x 2 ^- 2 + (i-l)(3j + l)x 2 ^- 3 + ^(j-l)(27j 3 -23j 2 -65j-6)x 2j - 4 + ... (2.5) 



The expression ( |2.4j ) is still valid for j = 0, but the factor is not polynomial, it reads 

P- 2 (x) = — — ^ — r = x~ 2 - x~ 3 + x~ 4 - x" 5 + (2.6) 
x(x + 1) 



and agrees with the large x expansion (2J3) for j = 0. 



This structure is strongly suggestive of algebraic recursion relations between the M's. 
We indeed found recursion relations between M^ 1 and the numbers of lower j, for the 



3 We have actually proved the cases j = 0, 1, 2, 3 recursively. Although some of them are quite 
tedious, all the proofs are based on the same strategy, which we will explicitly show for the case 
j = in the next section. The same strategy should extend to higher values of j. 
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first few j's. But the growing complexity of these relations when j is increased is not 
encouraging (typically, the structure of the recursion itself involves the detailed structure 
of the Mj meanders of order j). 

Semi— meanders. Like in the meander case, the number k of connected components of 
order n semi-meanders cannot exceed n. There is actually only one semi-meander of 
maximal number of connected components: it corresponds to having n concentric circles, 
each crossing the half-river exactly once, hence 

AfW = 1 . (2.7) 

The first values of MA ' are displayed in Table II below, for 1 < k < n < 14 (some of these 
numbers can be found in |J and H). 
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Table II: The numbers Mn of inequivalent semi-meanders of order n with k connected 
components, for 1 < k < n < 14, obtained by exact enumeration on the computer. 
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Like in the meander case, it is straightforward (a little less tedious here) to compute 
the first few M^~ J \ j = 0, 1, 2, ... and to guess their general structure as functions of n 
and j. We find 

M^ n) = 1 
M^" 1 ) = n - 1 

M(- 2 ) = l(n 2 + n-8) 

M^- 3 ) = l( n 3 + 6n 2 - 31n - 24) 



+ 25. 



n,4 



M^" 4 ) = l(n 4 + 14n 3 - 49n 2 - 254n) 

+ 155 n ,5 + 55 n , 6 (2.8) 
M^- 5 ) = _( n 5 + 25n 4 - 15n 3 - 1105n 2 - 1066n + 1680) 

+ 875 n , 6 + 425 nJ + 145 n , 8 
M^"- 6 ) = l( n 6 + 39n 5 + 145n 4 - 2895n 3 - 10226n 2 + 8616n + 31680) 

+ 4565 n>7 + 2925 n)8 + 1265 n , 9 + 425 n , 10 

M n - 7) = — in 7 + 56n 6 + 532n 5 - 5110n 4 - 50141n 3 - 20146n 2 + 377208n 
7! 

+ 282240) + 22345 n , 8 + 1720J n , 9 + 10085 n ,i + 396<5 n ,n + 1325 n , 12 . 

Note that in addition to some polynomial structure as functions of n, the numbers 
also incorporate some constant corrections for the j — 1 first values of n = j '• + 
l,j + 2, 2j — 1. The above suggests the following general form for the polynomial part 

of M^~ j) 



-Art + lj(3j - 5K'" 1 + - l)(27j 2 - 163j + 122)n^- 2 



pol J } - 2 24 (2.9) 

+ ^J(j - 1)(J " 2)(27j 3 - 354j 2 + 1163j - 1224)n 3 '" 3 + •••)• 



We also found the beginning of a pattern for the corrections to the polynomial part (^T 
of the form 



Cj5 n , 2 j-2 + 3cj5 n , 2 j-3 H ( 2 -10) 



where Cj = ( 2j ^ 1 ) / (2j + 1) is the Catalan number of order j, and the sum goes down to a 
6 n>j+ i term. 



Asymptotics. The data of the Tables I and II (see also [JTJ] for M n ) enable one to evaluate 
numerically the asymptotic behaviour of M n and M n which read respectively 

(12.25) n 
M n ~ const =775 — 

, (2-11) 

M n ~ const - — 4— 
n 2 

The exponent 7/2 for meanders was conjectured to be exact in 
Remarkably, the entropy of meanders 

s = lim -LogM n = Log R ~ 2.50 (2.12) 

is likely to be exactly twice that of semi-meanders 

8 = lim -LogM n = Log R ~ 1.25 . (2.13) 

n^oo 77, 

Considering M n (resp. M n ) as the number of foldings of a closed (resp. open with a 
fixed end) strip of 2n (resp. n — 1) stamps, general statistical mechanical considerations 
suggest that the leading behaviour of these numbers is exactly the same (i.e. oc R ), when 
expressed in the total number of stamps N = 2n (resp. N = n — 1) for large N. This leads 
to the relation s = # (= Log -R) between the corresponding ther mo dynamical entropies, 
i.e. R = R 2 between the leading terms. 



3. Arch statistics 




Fig. 5: A generic meander is the superimposition of two arch configurations. 

The most general meander of order n with arbitrary number of connected components 
is specified uniquely by its upper half (above the river) and lower half (below the river), 
as shown in Fig. 5. Both halves form systems of n non-intersecting arches connecting 2n 
bridges by pairs. Any two arches are either disjoint or included, one into the other. Any 
such system of arches will be referred to as an arch configuration of order n. 
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1 2j+2 



Fig. 6: Recursion principle for arch configurations. One sums over all 
positions 2j + 2 of the right bridge of the leftmost arch, which separates 
the initial configuration into two configurations X of order j and Y of order 
n — j — 1, respectively below the leftmost arch and to its right. 

3.1. Catalan numbers 

Let us first compute by recursion the number c n of arch configurations of order n, 
linking n pairs of bridges. Starting from one such arch configuration, let us follow the arch 
linking the leftmost bridge (position 1) to another, say in position 2j + 2 (see Fig. 6: this 
position has to be even). This arch separates the configuration into two sub-configurations 
of arches, X of order j (below the leftmost arch), and Y of order n — j — 1 (to the right 
of the leftmost arch). Summing over the position of the right-hand bridge of the leftmost 
arch (2j + 2), we get the following simple recursion relation 

n-l 

c n = Cn-x-j , n > 1 , (3.1) 

3=0 

where we have set cq = 1 . This rather simple example is nevertheless archetypical of a gen- 
eral type of reasoning used to count any relevant number associated to arch configurations. 
The scheme of Fig. 6 is quite general. 

The relation ( |3.1| ) is the defining recursion relation of the celebrated Catalan numbers 
c n , which count, among other things, the numbers of parenthesings (with n opening and 
n closing parentheses) of words of n + 1 letters. It follows immediately that 

2n + l \ n J ~ n!(n+l)! " K ' 

In the following, we will need the generating function C(x) of Catalan numbers, 



C(x) = £ 



71 = 

1 + x + 2x 2 + 5x 3 + Ux 4 + 42x 5 + 132x 6 + ■ ■ 



(3.3) 



subject to the algebraic relation 

xC(x) 2 = C(x) - 1 , (3.4) 

due to (|3.1|). One gets 

C(x) = . (3.5) 
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Fig. 7: A particular 5-component meander of order 5, and the corresponding 
arch configuration. 

3.2. Arch configurations and meanders 

Our interest in arch configurations is motivated by the fact that both meanders and 
semi-meanders can be built out of them. For instance, we get a direct one to one correspon- 
dence between n-component meanders of order n and arch configurations of order n: any 
arch configuration is completed by reflection wrt the river (see Fig. 7). As a consequence, 
we have 

= c n , (3.6) 



the total number of arch configurations of order n, in agreement with (|2.2| ). 

More generally, any multi-component meander of order n is obtained by superimpos- 
ing any two arch configurations of order n, one above the river, one below the river, and 
connecting them through the 2n bridges. As a consequence, we find the sum rule 



n 

£ 

k=i 



= (c n ) 2 , (3.7) 



expressing the total number of multi-component meanders of order n as the total number 
of couples of (top and bottom) arch configurations of order n. This is readily checked on 
the data of Table I: the sums of numbers by columns are equal to the square of the lowest 
number (Mn = c n ) in each column. 





(a) (b) 

Fig. 8: Any semi-meander (b) is obtained from the superimposition of 
an arbitrary arch configuration (top of (a)) and a rainbow arch configuration 
(bottom of (a)) connected through the In bridges (n = 5 here), and by folding 
the river as indicated, thus identifying the bridges by pairs. The process (a) 
<-» (b) is clearly reversible. 
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It is instructive to note that any multi-component semi-meander of order n may also 
be obtained as the superimposition of one arbitrary arch configuration of order n above an 
infinite river, and a "rainbow" arch configuration below it (see Fig. 8). The rainbow arch 
configuration of order n, denoted by lZ n , consists of n arches linking the opposite pairs 
of bridges: (1, 2n), (2, 2n — 1), ...,(n, n + 1). Note that the number of bridges is doubled 
in this representation. To recover the semi-meander, one simply has to fold the infinite 
river into a semi-infinite one, identifying the 2n bridges by pairs according to the rainbow 
arches as indicated in Fig. 8. As a consequence, we get the semi-meander version of the 
sum rule ( |3.7|) for meanders 

n 

£ = n , (3.8) 

k=l 

expressing the total number of multi-component semi-meanders of order n as the total 
number c n of arch configurations of order n (completed by the lower rainbow TZ n to yield 
the semi-meanders). This sum rule is readily checked on the data of Table II: the sums of 
numbers by columns are the Catalan numbers. 

The main difficulty in the meander and semi-meander problems is to find a direct 
way, for given arch configurations, to count the number of connected components of the 
resulting meander or semi-meander. This however is far beyond reach. Nevertheless, it 
is instructive to gather more refined statistical informations on the distribution of arches 
in random configurations, in view of a tentative generalization to the arch statistics of 
meanders with fixed number of connected components. 

3. 3. Statistics of arches 

Exterior arches. Let us first compute the distribution law of exterior arches. By 
exterior arches, we mean arches which have no other arch above them. For instance, the 
arch configuration of Fig. 7 has two exterior arches, a rainbow lZ n has only one exterior 
arch, etc... Let E(n : k) denote the number of arch configurations of order n with exactly 
k exterior arches. A simple recursion relation can be obtained in the same spirit as for 
Catalan numbers, by following the general scheme of Fig. 6. Starting from an order n 
arch configuration with k exterior arches, let us consider the leftmost arch, starting at the 
leftmost bridge. It is clearly an exterior arch. Let 2j + 2 be the position of the right bridge 
of this arch. Again, the arch separates the configuration into two arch configurations. The 
one below the arch is an arbitrary configuration among the Cj arch configurations of order 
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j. The one to the right of the arch is an arch configuration of order n — j — 1, with k — 1 
exterior arches. This leads to the recursion!! 



n-l 



E(n, k) = c j E ( n -3-l,k-l), (3.9) 

3=0 

for k > 1 and n > 1, with the initial condition E(n, 0) = an d E(n, k) = for k > n. 
This defines the numbers E(n,k) uniquely, and it is easy to prove that 

„, k {2n-k\ k(2n-k-l)\ 

E(n,k) = r = ttt^ ■ 3.10 

K ' J 2n-k\ n J n\{n - k)\ v ; 

Another way of characterizing the distribution of arch configurations according to 
their number of exterior arches is through its factorial moments, defined as 

where (J) = k(k - l)...(k - I + Thanks to the identity 

A/fc\^ n 21 + 1 f2n + l\ f2n\ ( 2n \ 

k=l 



which is proved in appendix A (note that for / = ( |3.12| ) gives the total number c n of arch 
configurations of order n), we deduce that 

, , , n\(n + l)\ 
\lP- = (a + 1) (n-0l(n + i + l)! - <3 ' 13> 

In particular, for 1 = 1, the average number of exterior arches for arch configurations of 
order n reads 

Wext = —~ ■ (3-14) 

n + 2 

In the limit of infinite order, n oo (the thermodynamic limit), the Zth factorial 
moment of the distribution of exterior arches tends to 21 + 1. In particular, the average 
number of exterior arches tends to 3. 

Interior arches. An interior arch of an arch configuration is an arch containing no 
other arch below it. For instance the configuration of Fig. 7 has 3 interior arches, a rainbow 
lZ n has only one interior arch, etc... Let I(n, k) denote the number of arch configurations 



One can also show that E(n, k) = E(n, k + 1) + E(n — 1, k — 1). 
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of order n with exactly k interior arches. As before (see Fig. 6), we can derive a recursion 
relation by following the leftmost exterior arch of any order n arch configuration with k 
interior arches. Three situations may occur: 

(i) the arch links the first bridge to the last one. Below this arch, we may have any 
arch configuration of order n — 1 with k interior arches. 

(ii) the arch links the first bridge to the second, hence it is an interior arch. To its 
right, we have any arch configuration of order n — 1, with k — 1 interior arches. 

(iii) the arch links the first bridge to the one in position 2j + 2, j = 1, — 2. It 
separates the configuration into two configurations of respective orders j (below the arch) 
and n — j — 1 (to the right of the arch), with respectively q and k — q interior arches, 
q = l,...,k- 1. 

This leads to the recursion relation 

n-2 fc-1 

I(n, k) = I(n -l,k) + I(n-l,k-l) + J2Yl r & ^ n -J-^k-q), (3.15) 

j=l q=l 

for n > k > 1, and with the initial condition 7(1, k) = 5k,i, and q) = for q > j. This 
allows one to solve for 

T( ,n _ ^( n \( n \ n\(n-l)\ 

[U ' K) ~ n\k)\k-l) ~ k\(k-l)\(n-k)\(n-k + l)\ ■ [ J 

With these numbers, we compute 



k=i 



A proof of this identity is given in appendix A. Hence the factorial moments of the distri- 
bution of interior arches read 

( W ,mt - EL, /(».*) 

(n\ ( 2n-l \ {O.IO) 
VU \n-l + l) 

(„ 2 ") ' 

For instance the average number of interior arches on an arch configuration of order n is 

14 



Fig. 9: The S duality between arch configurations is indicated by arrows, 
for orders n = 1, 2, 3. 



From (|3.18| ), it is easy to see that the Zth factorial moment behaves as 

<©>" ~ M <X20) 

for large n. All these numbers tend to infinity when n is large, hence diverge in the 
thermodynamic limit. 

A last remark is in order. The numbers J(n, k) (|3.16 ) satisfy a symmetry relation 



I(n,k) = I(n,n+l-k) , (3.21) 

which is in fact the expression of a duality between arch configurations. Singling out the 
leftmost exterior arch, any arch configuration may be expressed as on Fig. 6. The duality 
transformation is defined recursively as 

5(0) = 

^ ^ (3-22) 

S((X\Y) = fsm\ S(X) 

This transformation is involutive (505 = 1) and exchanges disjoint and nested arches 
as exemplified in Fig. 9. The transformation S changes the number of interior arches of a 
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given configuration of order n from k to n + 1 — k, as is readily shown by recursion. This 



explains the symmetry (|3 . 2 1 



Arches which are both interior and exterior. Let us compute the total number 
EI (n) of arches which are both interior and exterior in all the arch configurations of order 
n. We have EI(1) = 1. By the recursion process of Fig. 6, we can relate the number 
EI {n) to the numbers EI(j), j < n in the following manner. Summing over the position 
2j + 2 of the right bridge of the leftmost exterior arch in any order n configuration, we 
find that only the configuration Y of order n — 1 — j, to the right of this arch, contributes 
to EI {n) through the total of EI(n — 1 — j) arches which are both interior and exterior. 
This number comes cj times, corresponding to the arbitrary choice of the configuration X 
of order j below the leftmost exterior arch. When j = 0, the leftmost arch is both interior 
and exterior, and there are c n _i such configurations, hence in addition to the EI(n — 1) 
arches which are both interior and exterior, we also have c n _i leftmost arches, which are 
both interior and exterior. There is no contribution from j = (n — 1) (for n > 2) as the 
only exterior arch cannot be interior. This yields 

n-2 

EI {n) = EI(n- j - 1) + EI(n - 1) + c n _i , (3.23) 

i=i 

valid for n > 2. If we define EI(0) = 1, this takes the form 

71-1 

EI(n) = J2 c j EIin-1-j) , (3.24) 

j=0 

hence 

EI(n) = c n . (3.25) 

For all n, there is therefore exactly one arch in average which is both exterior and interior. 

Arches of given depth. The depth of an arch in an arch configuration is defined as 
follows. The exterior arches are the arches of depth 1. The arches of depth k + 1 are the 
exterior arches of the arch configuration obtained by erasing the arches of depth 1, 2, ...,k. 
For instance the rainbow configuration 1Z n of order n has one arch of each depth between 
1 and n; the arch configuration of Fig. 7 has 2 arches of depth 1, and 3 arches of depth 2, 
which are also its interior arches. Let us compute the total number A(n, k) of arches of 
depth k in all the arch configurations of order n, by the usual reasoning of Fig. 6. Starting 
from any order n arch configuration, consider as usual the leftmost arch. It links the first 
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bridge to the one in position 2j + 2, j = 0, n — 1. The arch separates the configuration 
into one of order j (below the arch, X in Fig. 6), whose depth k — 1 arches give rise to 
depth k arches of the initial configuration, and one of order n — j — 1 (to the right of the 
arch, Y in Fig. 6), whose depth k arches are depth k arches of the initial configuration. 
This translates into the recursion relation 



n-l 



A(n, k) = }^ (A(j, k - 1)^-1 + A(n k) Cj ) , (3.26) 

for n > and k > 2, with A(n, k) = if k > n. The initial condition is the total number 
A(n 1 1) of external arches of all order n arch configurations^, namely 

A(„,l) = ±mn,i) = ^f 2 ;^ 1 ) , (3.27) 

j=l ' 

according to (|3.12| ) for I = 1. This fixes the A(n, k) completely, and we find 

,/ ,n 2k + lf2n + l\ ( 2n \ ( 2n \ /o ooN 

This is nothing but the rhs of the identity ( |3.12j ), with the substitution I — ► k (indeed 
we have taken the I = 1 value of this expression as initial condition). We still lack of a 
good explanation for this phenomenon. Note that the formula ( |3.28| ) can be used as the 
definition of A(n, 0) = c n , the total number of arch configurations, thus interpreting the 
depth arch number to be 1 for each configuration. We can define the average number of 
depth k arches as 

These numbers coincide with the factorial moments (|3.13|) of the distribution of external 
arches. In particular, in the thermodynamic limit (n — > oo), the average number of depth 
k arches tends to 2k + 1. So the average arch configuration has 3 exterior arches, 5 arches 
of depth 2, 7 arches of depth 3, etc ... 

The total number of arches on a given configuration is n. Therefore, we have the 
following sum rule 

oo 

J2 A (n,j) = nA(n,0) . (3.30) 



Indeed the initial condition A(l, 1) = 1 would be sufficient, as the expression is recur- 

sive wrt k + n. 
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The above result for the thermodynamic averages enables us to define an average maximal 
depth in arch configurations of order n, denoted by d(n), by 

r Mn,3) 
n = > lim — 

^ n^oo A(n, 0) 

J ~ (3.31) 

a(n) 

= J2 ( 2 ^' + 1 ) = d(n)(d(n) + 2) . 

This gives 

d(ra) = v / n+T- 1 . (3.32) 



4. Recursion relations for semi— meanders 

In the following, we will mainly concentrate our efforts on semi-meanders. However, 
we will mention whenever possible the extensions of our results to meanders. 

4-1. The main recursion 

Let us present a simple algorithm for enumerating the semi-meanders of order n with 
k connected components. In this section, we will work in the infinite river /lower rainbow 
arch-framework for semi-meanders (see Fig. 8 (a)), namely consider a semi-meander as the 
superimposition of the lower rainbow configuration lZ n , and some upper arch configuration. 

Starting from a semi-meander of order n with k connected components, let us con- 
struct several semi-meanders of order n + 1 with k connected components in the following 
way, as indicated in Fig. 10 (I) (i)-(ii). 

(i) Pick any exterior arch of the upper arch configuration of the semi-meander and 
cut it. 

(ii) Pull the two edges of the cut across the river (the left part of the exterior arch 
to the left, the right one to the right), and paste them around the lower rainbow, thus 
increasing the rainbow configuration by one arch (1Z n — > 1Z n +i), and the number of bridges 
by 2. 

The result is a certain semi-meander of order n + 1 with the same number k of 
connected components. So for each semi-meander M. of order n with k components, 
we can construct E{M) semi-meanders of order n + 1 with k components, where E{M) 
denotes the number of exterior arches of M. All these new semi-meanders are clearly 
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Fig. 10: The construction of all the semi-meanders of order n + 1 with arbi- 
trary number of connected components from those of order n. Process (I): (i) 
pick any exterior arch and cut it (ii) draw its edges around the semi-meander 
and paste them below. The lower part becomes the rainbow configuration 
1Z n +i of order n + 1. The process (I) preserves the number of connected 
components. Process (II): draw a circle around the semi-meander of order n. 
The process (II) adds one connected component. 

distinct. There is however another way of generating more semi-meanders of order n + 1, 
indicated in Fig. 10 (II). Starting with any semi-meander of order n with k — 1 connected 
components, one just adds an extra circular loop around it, which transforms the lower 
rainbow of order n into 7Z n +i, and adds 2 bridges. The resulting semi-meander of order 
n + 1 has clearly k connected components. Such a semi-meander cannot be obtained from 
some order n semi-meander by the procedure (I), because it has only one exterior upper 
arch, whereas the process (I) produces at least two exterior arches. 

So all the order n + 1 semi-meanders constructed by (I)- (II) are distinct. Conversely, 
given a semi-meander of order n + 1 with k components, two cases may occur: 

(a) it has only one exterior upper arch. In this case, it is surrounded by one circle, and 
therefore arises from the order n semi-meander with k — 1 components inside the circle, 
through (II). 

(b) it has at least two exterior upper arches. Cutting the lower exterior arch of the 
lower rainbow TZ n+ i, pulling the edges of this arch above the upper configuration (the left 
edge by the left side, the right one by the right side, both suppressing the left and right- 
most bridges) , and finally pasting the two edges on the upper side of the river (thus creating 
an exterior arch on the upper configuration), one gets an order n semi-meander with k 
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components, which leads to the initial semi-meander of order n + 1 through (I)-(i)-(ii), 
using the exterior arch constructed above. 

This proves that the procedures (I)- (II) give a recursive algorithm for constructing all 
the semi-meanders of order n + 1 with k connected components from the k and k — 1- 
component meanders of order n. More precisely, these procedures can also be viewed as 
a recursive algorithm for constructing all the arch configurations of order n from those of 
order n — 1. When completed by a lower rainbow configuration lZ n , so as to give multi- 
component semi-meanders, we have the interesting property that (I) preserves the number 
of connected components of the semi-meander, while (II) increases it by 1, which allows 
us in principle to follow this number of components throughout the construction. 

The case k = 1 of connected semi-meanders is special, in the sense that the corre- 
sponding arch configurations are obtained by successive actions of the procedure (I) only. 
As a consequence, the number of connected semi-meanders of order n + 1 is equal to the 
total number of exterior upper arches of all the semi-meanders of order n. 

4-2. Expression in terms of arch numbers 

We now want to translate into some recursion relation the transformations (I)— (II) . 
It is best expressed in terms of the numbers of arches of given depth in the upper con- 
figurations of the semi-meanders. Recall that to build a semi-meander, one just has to 
close an upper arch configuration using a lower rainbow arch configuration. Therefore, 
we can decompose the set of arch configurations of order n into a partition, according to 
the number k = 1,2, of connected components of the corresponding semi-meander. 
Denoting respectively by A n and SM n ^ the sets of arch configurations of order n and of 
semi-meanders of order n with k connected components, we have therefore 

An = UJU SM n , k . (4.1) 

The statistics of arches extensively studied in sect. 3 only give some global information on 
semi-meanders, irrespectively of their number of connected components. Taking the latter 
into account is a subtle refining of the above study. In particular we do not expect the 
statistics of arches to be the same when the number of connected components is fixed. We 
rather want to use sect. 3 as a guideline for the study of semi-meander arch statistics. In 
particular, we will use analogous definitions for observable quantities, and derive recursion 
relations between them. 
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For a given semi-meander M. of order n with k components, let us denote by A(A4,j) 
the number of arches of depth j in its upper configuration, and let 

A k (n,j) = £ A(M,j) , (4.2) 

MeSM n , k 

denote the total number of arches of depth j of the upper configurations of the semi- 
meanders of order n with k connected components. To make the contact with the numbers 
A(n, j) of arches of depth j in the arch configurations of order n (eq. (|3.28"|) ), let us mention 
the obvious sum rule 

n 

MnJ) = A(nJ) . (4.3) 

k=i 

This simply states that summing the numbers of upper arches of depth j over all (multi- 
component) semi-meanders of order n amounts to counting the total number of arches 
of depth j in all the arch configurations of order n. This sum rule is parallel to the set 



decomposition (|4.1|) , and simply evaluates the number of depth j arches on both sides. 

To write a recursion relation, we must follow the evolution of the depth of arches 
through the processes (I) and (II). In the process (II), it is clear that the arches of depth 
j — 1 in the order n semi-meander become arches of depth j in the order n + 1 one, hence 
we get the contribution 

A k _ 1 (n,j-1) = Yl A (MJ-1) (4.4) 

-M6SM re , fe _i 

to the total number of arches of depth j. The process (I) of Fig. 10 is more subtle. We start 
from a meander of order n with k components, Ai n ,k- For each of its E(Ai n ,k) exterior 
arches, the process (I) affects the depth of its arches in the following way. The arches 
of depth j + 1 below the specified exterior arch (configuration denoted by Y on Fig. 10) 
become arches of depth j. The arches of depth j — 1 to the left (X) and to the right (Z) 
of this exterior arch become arches of depth j. Summing over all the A(A4 nj k, 1) exterior 
arches of A4 n ,k, we get the contribution 

A(M n , k J + 1) + [A(M n , k , 1) - 1] A(M n , k , j - 1) (4.5) 



to the total number of arches of depth j. Adding the contributions (|4.4|) and ( fl.5|) of the 



processes (II) and (I), we finally get the main recursion relation, for k > 2 and j >2 
A k (n + l,j) = A k - 1 (n,j-l)+A k (n,j + l)+ ^ [A(M, 1) - 1] A(M,j - 1) . (4.6) 

MeSM n u 
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Like in the arch configuration case, we can introduce the number of arches of depth in a 
semi-meander, in such a way that 

A(M,0) = 1 VMe SM n , k . (4.7) 

In particular, 

A k (n,0) = card(SM n , fc ) = . (4.8) 
For j = 1, the recursion reads 

A k (n + 1,1) = A fc _i(n, 0) + A k (n, 2) + 2A k (n, 1) . (4.9) 

For j = 0, the recursion relation becomes 

A fc (n + 1,0) = A fc _i(Ti,0)+A fc (n,l) . (4.10) 



Eqs.(0Q(g are all valid for k > 2. 

The case fc = 1 is different since it uses the procedure (I) only. However, the above 
equations remain valid upon setting Ao(n, j) = for any n, j. This leads to 

A^n + lJ) = A 1 {n,j + 1)+ ]T [A(M, 1) - 1] A(M, j - 1) , (4.11) 

M6SM nil 

valid for j >2. For j = 1,0, we have 

A l (n + 1,1) = A 1 (n,2) + 2A 1 (n,l) 

(4.12) 

Ai(n + 1,0) = Ai(n,l) . 

The last equation expresses the fact that the total number M n +i of connected semi- 
meanders of order n + 1 is equal to the total number of exterior (depth 1) arches on all 
the semi-meanders of order n. 

The abovementioned leading behaviour ( |2.11|) of the connected semi-meander number 

M n oc R n , (4.13) 

enables one to recover the value of R from the thermodynamic limit of various averages 
over semi-meanders. 
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We now denote by 

the average of a function / over the semi — meanders of order n with k connected compo- 
nents. Dividing eq.(fOI) by Ak(n, 0), we find a relation between averages 



A k (n+ 1,0) A fc _i(n,0) . 

— — — — — (A(M,j)) n +i,k = , ( ' (A(M,j -l)) n ,fc-i 
A k (n,0) A k (n,0) (4.15) 

+ (A(M,j + 1) + [A(M) - l]A(M,j - l)) n>fc , 
for k > 2. For /c = 1, we have 

^r^ mjiw = mi+ij+^-wi-i)),,!, (4.i6) 

Ai[n, 0) 
for j > 2 and for j = 1 , 

Al ( n ( +1 ' 0) (A(M,l))n+i,i = (A(M, 2) + 2A(M, l)) n ,i 

Ai (ra + 1,0) . . . . . iN . v ; 

i,K0) = WM 1)> - 1 ' 

This last equation enables us to interpret R as the average thermodynamic number of 
exterior arches 

R = lim + = Urn %±i = HWDW . (4.18) 

^.5. Mean field approximation 

In the thermodynamic limit, in view of the result for the statistics of arches, we 
expect all these averages to tend to some finite values. A standard approximation when 
dealing with averages over some statistical distribution is the mean field approximation. 
It expresses that the observable quantities (here the arch numbers) are fixed to their mean 
value. In particular, in this approximation, the term correlating the arches of depth 1 and 
j — 1 is replaced by 



([A(M) - l]A(M,j - l)) n>k — ([A(M) - l]) nik {A(M,j - l)) n>k . (4.19) 
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From now on, we focus our attention to the connected case k = 1, where this approximation 
gives the system 

Al{ " L ° ) (A(M,j)) n+1 , 1 = (A(M,j + l)) n ,i + ([A(M) - l]) ntl {A(M,j - l)) njl 



A t (n, 0) 
A 1 (n + 1,0) 
Ai(n,0) 

A x (n + 1,0) 



(A(A*,l)>n+i,i = (A(M,2) + 2A(M,l))n,i 

(A(M,l))n,l ■ 



Ai(n,0) 

(4.20) 

Starting from the only semi-meander of order n = 1, (a circle crossing the river through 
two bridges), for which 

{A{M,j))x, x = 6 jA , (4.21) 
we use the system (|4.20[ ) recursively to get the quantities 

a n (j) = (A(M,j)) n ,i , (4.22) 

for any n > 1. Let us prove that when n — > oo they converge to the values 

(A(M,j))oo,i = 3 Vj . (4.23) 



If the quantities ( [4.22| ) converge, their limits, denoted by a(j) = {A(JA, i))oo,ij are 



subject to the recursion relation (as fixed points of the system ( |4. 2U| ) ) 

R a(j) = a(j + 1) + [a(l) - l]a(j - 1) , (4.24) 



valid for j > 2. On the other hand, a(l) = i?, from eq.( |4.18| ). We are thus left with the 
recursion relation 

a(j + 1) = £a(i) - (fl - l)a(j - 1) , (4.25) 

for j > 2, whereas 

a(2) = a(l) (i?- 2) . (4.26) 

Note that necessarily R > 2 for the average number a(2) to be positive. With the initial 
conditions a(0) = 1 and a(l) = R, we find in the 1-component case (k = 1) that 

a(j) = ia + ^-Sp-r 1 )- (4-27) 
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This is the mean field solution for the numbers of depth j arches in connected semi- 
meanders. These numbers are further constrained by the fact that the average numbers of 
arches have to be non-negative, hence 



R > 3 . (4.28) 

Let us now prove that R < 3. For this purpose, let us show by recursion that the 
inequalities 

3 > (A(M, l))n,i > (A{M, 2)) n ,i > • • • (4.29) 

are satisfied for all n. The inequalities ( |4.29|) are obvious for n = 1. Suppose (|4.29|) is true 
for some value of n, then using the recursion relations ( |4.20[ ), and performing the difference 
of any two recursion relations for consecutive values of j, we get 

o n (l)o n+ i(l) < 3 a n (l) 
a n (l)(o n+ i(2) - a n+1 (l)) < a n (l)(o n (l) - 3) (4.30) 
a n (l)(a n+1 (j + 1) - a n+1 (j)) < (a„(l) - l)(a n (J) - a n (j - 1)) 

which imply ( J4.29[ ) at order n + 1. In particular we always have a n (l) < 3, hence in the 
large n limit, a(l) = R < 3. 

This leaves us with the unique possibility 

R = 3 => a(j) =3 Vj > , (4.31) 

which is the actual limit of the mean field approximation to arch numbers of given depth. 
Note that in this case the recursion Q4.27 ) is satisfied for j = 0, 1 as well. Of course this is 



only a rough approximation of the exact R observed in ( |2.11| ), as is usually expected from 
a mean field approximation. 

4-4- Mean field is wrong for arch configurations 

The recursion relation (fO|) can be used to recover some information on arch con- 
figurations, by summing it over all the numbers k of connected components. Recall that 
J2k Ak(n,j) = A(n,j), the total number of arches of depth j in the arch configurations of 
order n. Hence, summing over k = 1, 2, n + 1 in ( |4.6| ), we get 

A(n + l,j) = A(n,j-1) + A(n,j + 1)+ ^ [A(M, 1) - 1] A(M,j - 1) , (4.32) 

Me An 
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where A n denotes the set of arch configurations of order n. Using the explicit formula 
(|3.28|) for the number of arches, we deduce that 



[A(M,l)-l}A(M,j-l) = 2A(n,j). (4.33) 

MeA n 

This is obtained by noticing that the numbers A{n, j) satisfy the linear recursion relation 
A(n + = A(n, j - 1) + A(n, j + 1) + 2A(n, j) , (4.34) 
for n, j > 1, as a consequence of the expression (|3.28|) of A(n, j) as a difference between two 



binomial coefficients. This gives the value of the correlation function between the numbers 
of arches of depth 1 and those of depth j — 1 in all arch configurations. Dividing ( |4.33j ) 
by the number of arch configurations A(n, 0) = c n , we get a relation for averages. If we 
denote by 

< /(A1)> " = IRo) £ }(M) • (4S5) 

v ' MeA n 

this relation reads 

([A(M, 1) - 1] A(MJ - l)) n = 2(A(MJ)) n . (4.36) 

In the thermodynamic limit n — > oo, as (A(A4, j)) n — > 2j + 1 (see eq.( |3.29|) ), we find that 
([A(M, 1) - 1] A(M,j - l)) n - ([A(M, 1) - l]) n (A(MJ - l)) n 

(4.37) 

— > 2(2j + l)-(3-l)(2j-l) = 4. 

n — >oo 

We see that the mean field approximation used in the previous section for semi-meanders 
is clearly wrong in the case of arch configurations. Note also that the result ( |4.37|) is 
independent of j. 

However, it is easy to solve the mean field equations for arch configurations, which 
lead to a growth (2 + V2) n for the total number of arch configurations, whereas the exact 
behaviour is 4 n (as readily seen from the large n expression of the Catalan number c n by 
use of the Stirling formula). 
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4-5. Improved mean field approximation for semi-meanders 

This suggests to try a deformed thermodynamic mean field-type ansatz for the num- 
bers of arches in semi-meanders, of the form 



([A(M) - l]A(M,j - l)) n ,k — K k + ([A(M) - l])^ k (A(M, 3 - 1)W , (4-38) 



n — >oo 



where K k are some constants to be determined. For k = 1, the main recursion relations 
Q4.11|) and ( |4. 12|) read, in the thermodynamic limit 



R ai (j) = ai (j + 1) + K! + [ ai (l) - l] ai (j - 1) fori>2, 
ai(2) = ai(l)GR-2) 

(4.39) 

ai(l) = R ai(0) 
ai(0) = 1 . 

We determine the constant K\ by imposing that the main recursion be satisfied for the 
value j = 1 as well, in which case we get 

R ai(l) = R 2 = R(R-2) + K x + (R - 1) => K 1 = J R+1. (4.40) 

Next the recursion is easily solved by noticing that 

ai (j + 1) - a(j + 1) = R Mi) - o-j) - (£ - 1) ( ai (j - 1) - a(j - 1)) , (4.41) 

where 

a=|±l. (4.42) 



This is exactly the same recursion as in the mean field case ( [4. 25[ ) , except that in the new 
variables ot\(j) = a\(j) — aj, the initial conditions read 

ai(0) = l ai(l)= ^_ 2 ■ (4.43) 



The solution reads 



R + l . 1 



«i(i) = + ( ^^ 2)2 (3 + (^ 2 -4i?+l)(i?-iy) . (4.44) 

Like in the mean field case, the solution should not grow exponentially with the depth j. 
This imposes that the prefactor of (R — 1) J in ( f4.44j ) vanishes, i.e. 

R 2 -4R + l = # = 2 + ^ = 3.732... (4.45) 
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hence 

ai (j) = (1 + V3)j + 1 . (4.46) 

The value of R found in Q4.45 ) is now above the observed value (|2.11|) . 

When applied to the arch configurations, this improved mean field ansatz turns out 
to yield the exact thermodynamic solution 



(A(M,j)) c 



2j+l 



(4.47) 



which lead to the correct growth 4 n for the number of arch configurations of order n. 

It is possible that even more refined improvements of the mean field ansatz give better 
approximations for semi-meanders. One can for instance think of replacing the constant 
K\ above ( [4.38| ) by some specific function of j, but such a function is no longer fixed by 
initial conditions. 



4-6. Other approximations using the main recursion 

Starting from a particular arch configuration, let us generate its "descendents" by 
repeatedly using the processes (I) and (II) of Fig. 10. After a given sequence of n ((I) or 
(II)) steps, let us denote by (to, e) respectively the total number of arch configurations 
generated and their average number of exterior arches. Then let us adopt the following 
"mean field type" recursive algorithm for an extra action by (I) or (II) 

(to, e) — > (me, e) 

{I) (4.48) 
(™>e) — > (m,l) . 

The second line of ( |4.48|j is clear, as (II) builds one arch configuration out of each inital arch 
configuration, with only one exterior arch. The first line incorporates two suppositions. 
First, the number of arch configurations generated is approximated by a mean field value, 
namely the product of the initial number m of arch configurations by the average number 
e of exterior arches. Second, the average number of exterior arches is supposed to be 
unchanged. 

Let us now choose as starting point (to = 1, e) the upper arch configuration of a semi- 
meander. Repeated actions of (I) only will generate its descendents which are themselves 
semi-meanders. Implicitly, ( |4.48| ) supposes a certain number of properties on the starting 
point of the recursion, like for instance the fact that repeated actions of (I) keep the average 
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number of exterior arches e unchanged, while the number of semi-meanders obtained is 
multiplied by e each time. Thus in this scheme, e will be identified with R, governing the 
large n behaviour of the number of semi-meanders M n ~ R n . The scheme (|4.48| ) is only 
intended as an approximation and for that purpose the initial semi-meander is supposed 
to be very large, with a number e of exterior arches equal to the average R. Let us try 
to evaluate e by identifying the total number p n (e) of configurations generated after all 
possible sequences of n actions of (I) or (II) with the total number c n ~ 4 n (this assumes 
that c n gives the correct large n behaviour of the total number of descendents of a given 
configuration, independently of this configuration) . We start with only one semi-meander, 
hence po(e) = 1. After one step, according to eq.( |4.48j ), we have generated pi(e) = e + 1 
arch configurations. Let us decompose the number p n {e) into 

Pn(e) = P W(e) + pW(e) , (4.49) 

where we make the distinction between the total number of arch configurations obtained 
by the process (I) (resp. (II)) in the last step ( f4.48|) from the p n _i(e) previous ones. The 
algorithm (|4.48| ) leads to the recursion relations 

P { n\e) = ep^Ue) + pS(e) 
P ( n n) (e) = p&(e) + pS(e). 

In terms of the vectors P n (e) = (pn 7 ' ) (e),pi / ' f ' ) (e)) t , this takes the matrix form P n (e) = 
M(e)P n _i(e), with 

M(e) = [\ ]) . (4.51) 



x l 1 1 

For the total number of arch configurations generated to behave like 4 n , we simply have 
to write that the largest eigenvalue of this matrix is 4, namely that the determinant of 
M(e) — 41 vanishes (this also fixes the other eigenvalue of M(e) to be 2/3). This gives 

e = — = 3.666... (4.52) 

This estimate of R is the closest we can get to the numerical estimate ( |2.11| ). But the 
sequence of approximations used to get ( |4.52j ) should certainly be refined. 

A last remark is in order. One could wonder how much the above estimate depends 
on the starting point of the algorithm. In particular, starting with any semi-meander 
with a finite number k of connected components and a large order, we end up with the 
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same estimate ( |4.52|) for the average number of exterior arches. This in turn infers an 
estimate for Rk = Ri = R = e = ll/3 for the numbers Rk governing the large order 
behaviour of Mn ~ (Rk) n - Moreover the average number of exterior arches for all the 
arch configurations (generated in our scheme by both (I) and (II)) reads, after n steps, 

(I) , , (II) 

(ext) n = — ^ tjjt — . (4.53) 

Pn +Pn 

For large n, the vector (pn\pn^) tends, up to a global normalization, to the eigenvector 
(3, 1) associated to the largest eigenvalue A max = 4 of the matrix M(ll/3) ( fl51| ). Therefore 



the average number of exterior arches for all arch configurations is estimated as 

(ext) = = 3 , (4.54) 



which coincides with the exact value, as given by ( |3.14| ). 



5. Matrix model for meanders 

Field theory, as a computational method, involves expansions over graphs weighted by 
combinatorial factors. In this section, we present a particular field theory which precisely 
generates planar graphs with a direct meander interpretation. The planarity of theses 
graphs is an important requirement, which ensures that the arches of the meander do not 
intersect each other, when drawn on a planar surface. The topology of the graphs in field 
theoretical expansions is best taken into account in matrix models, where the size N of the 
matrices governs a topological expansion in which the term of order j\T 2-2?l corresponds 
to graphs with genus h. The planar graphs (with h = 0) are therefore obtained by taking 
the large iV limit of matrix models (see for instance H for a review on random matrices). 

5.1. The matrix model as combinatorial tool 

Random matrix models are useful combinatorial tools for the enumeration of (con- 
nected) graphs . Typically, one considers the following integral over Hermitian matrices 
of size N x N 

Z(g,N) = J 5 (5.1) 
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where the integration measure is 



N 

dM = Yl dM u J[ dReMij dlmM i:j . (5.2) 

The rules of Gaussian matrix integration are simple enough to provide us with a trick for 
computing the expression ( |5.1| ) as a formal series expansion in powers of g. For instance, 

JdM e- NT ^M l3 M kl 8 u S jk 



(MijMki) 



auss 



fdMe- 1 ™^ N 
Expanding Z(g, N) in powers of g, we are left with the computation of 



V=0 



V\ 



(5.3) 



(5.4) 




Fig. 11: A ribbon graph with V = 2 vertices, E = 4 edges, and L = 4 
oriented loops. 



The result (|5.3|) suggests to represent any term in the small g expansion as a sum over 
graphs constructed as follows. 



§i,l 5j,k/N 



'4 J 4 



'l' J 2 l 2' J 3 1 3' J 4 1 4 ,3 1 



Fig. 12: The Feynman rules for the one matrix model: the matrix indices 
are conserved along the oriented lines, which form closed loops. Edges receive 
a factor 1/N, vertices (Ng), and all the matrix indices have to be summed 
over. 
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The graphs considered are ribbon graphs, i.e. with edges made of double-lines, oriented 
with opposite orientations (see Fig. 11 for an example). These two oriented lines represent 
the circulation of matrix indices i, j = 1, 2, A. Namely each end of an edge is associated 
with a matrix element My, the index i being carried by the line pointing from (resp. j 
by the line pointing to) this end of the edge, as shown in Fig. 12. An edge is therefore 
interpreted as the propagator between the states sitting at its extremities, according to 
( |5.3|) , hence each edge will be weighed by a factor 1/A (hence an overall factor N~ E , where 
E is the total number of edges of the graph) . To compute the term of order V in the small 
g expansion ( |5.4j ), one simply has to sum over all the ribbon graphs connecting the V four- 
valent vertices corresponding to the V terms TrM 4 (an obvious 4-fold cyclic symmetry 
absorbs the factors 1/4). Each vertex has to be weighed by a factor Ng. These Feynman 
rules are summarized in Fig. 12. An overall weight also comes from the summation over all 
the matrix indices i = 1, 2, N running on the oriented loops of the graph. This gives a 
global factor N L for each graph, where L is the number of loops of the graph. For instance, 
the ribbon graph of Fig. 11 receives a total weight (Ng) 2 x A -4 x A 4 = (Ng) 2 . 

The second trick is the fact that this sum can restricted to connected graphs only, by 



taking the logarithm of the function ( |5.1|) 



F(g,N) = Log Z(g, A) = £ g v N v ~ E+L x ——— , (5.5) 

conn, graphs T 

where |Aut(r)| denotes the order of the automorphism group of T, i.e. the number of 
permutations of its (supposedly labelled) vertices leaving the graph invariant. This sym- 
metry factor results from the incomplete compensation of the factor 1/V\ by the number of 
equivalent graphs (with different labelling of the vertices). Finally, we identify the power 
of A as the Euler-Poincare characteristic of the graph Y 

X (Y) = V-E + L = 2-2/i, (5.6) 

which can be taken as the definition of the genus h of the graph. The number of oriented 
loops is indeed equal to that of faces F of the cellular complex induced by the graph, 
hence we can use the more standard definition of the Euler-Poincare characteristic x = 
V-E + F = V-E + L. 



Various techniques for direct computation of the integral (|5.1| ) have made it possible 
to enumerate connected graphs with arbitrary genus, and derive many of their properties. 
In this section, we consider a matrix model adapted to the meander enumeration problem. 
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Fig. 13: A sample black and white graph. The white loop is represented in 
thin dashed line. There are 10 intersections. 

5.2. The model 

The enumeration of (planar) meanders is very close to that of 4-valent (genus 0) 
graphs made of two self-avoiding loops (say one black and one white), intersecting each 
other at simple nodes The white loop stands for the river, closed at infinity. The black 
loop is the road. Such a graph will be called a black and white graph. An example is given 
in Fig. 13. The fact that the river becomes a loop replaces the order of the bridges by a 
cyclic order, and identifies the regions above the river and below it. Hence the number 
of meanders M n is 2 x 2n (2 for the up/down symmetry and 2n for the cyclic symmetry) 
times that of inequivalent black and white graphs with 2n intersections, weighed by the 
symmetry factor l/|Aut(r)|. The same connection holds between Mn and the black and 
white graphs where the black loop has k connected components. 






Fig. 14: A particular black and white graph with 6 intersections, and its 
two associated meanders. The automorphism group of the black and white 
graph is TL§. 
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For illustration, we display a particular black and white graph V in Fig. 14, together with 
its two corresponding meanders of order 3. The automorphism group of this black and 
white graph is Z6, with order Aut(r) = |Zg| = 6. The two meanders come with an overall 
factor 1/(2 x 6), hence contribute a total 2 x 1/12 = 1/6, which is precisely the desired 
symmetry factor. 

In analogy with the ordinary matrix model (|5.1| ), a simple way of generating black 
and white graphs is the use of the multi-matrix integral (with m + n hermitian matrices 
of size iV denoted by B and W) 

Z(m,n,c,N) = — / TT dB& FT dW^e~ N ^ p&»\w™) (5 . 7) 



a=l 0=1 



where the matrix potential reads 



p{bw,wW) = Y,^^+Y,^^-lY, B{(X)w{p)B{a)w{p) ' ( 5 - 8 ) 

a (3 a,/3 

and the normalization constant kn is such that Z(m,n : c = 0,N) = 1. In the following, 
the a and f3 indices will be referred to as color indices. 



a ; — < — . a' 

J K AY 

Sy Sj, k S^a'/N « j 1 < ',', < i 3 «' 



A* 



i > 1 „ • 
P j 1 k |S 

8i,i 8j, k Sp^p'/N 



•4p'4 

Nc8i i Sj j Si j S| : 8 a a >SR r' 

l' J 2 '2' J 3 '3' J 4 '4' J 1 ' ^'^ 



Fig. 15: The Feynman rules for the black and white matrix model. Solid 
(resp. dashed) double-lines correspond to black (resp. white) matrix ele- 
ments, whose indices run along the two oriented lines. An extra color index a 
(resp. 0) indicates the number of the matrix in its class, B^ a \ a = 1, 2, m 
(resp. j3 = 1,2, n). The only allowed vertices are 4-valent, and have 

alternating black and white edges: they describe simple intersections of the 
black and white loops. 



Like in the case ( |5.1|) , the logarithm of the function (|5.7| ) can be evaluated pertuba- 
tively as a series in powers of c. A term of order V in this expansion is readily evaluated 
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as a Gaussian multi-matrix integral. It can be obtained as a sum over 4-valent connected 
graphs, whose V vertices have to be connected by means of the two types of edges 

black edges {[B^^B^u) = ^^<W 

(5 9) 

white edges ([W^j [W^] kl ) = ^ V > 

which have to alternate around each vertex. The corresponding Feynman rules are sum- 
marized in Fig. 15. This is an exact realization of the desired connected black and white 
graphs, except that any number of loops! of each color is allowed. In fact, each graph 
receives a weight 

N 2 ~ 2h c v m b n w , (5.10) 

where b (resp. w) denote the total numbers of black (resp. white) loops. 

A simple trick to reduce the number of say white loops w to one is to send the number 
n of white matrices W to 0, and to retain only the contributions of order 1 in n. Hence 

f(m, c, AO = lim ^Log Z(m, n, c,N) = V N 2 ~ 2h c v m h ] . 

n^o n Autll ) 

b.Szw. conn, graphs Y \ \ / \ 

with one w loop 

(5.11) 

If we restrict this sum to the leading order A^ 2 , namely the genus contribution 
(h = 0), we finally get a relation to the meander numbers in the form 

fo(m,c) = lim —^f(m,c,N) 

00 n 2 P p (5-12) 

P =i F k=i 



where the abovementioned relation between the numbers of black and white graphs and 
multi-component meanders has been used to rewrite the expansion ( |5.11|) . 

5. 3. Meander numbers as Gaussian averages of words 

The particular form of the matrix potential (|5.8| ) allows one to perform the exact 
integration over say all the W matrices (the dependence of P on W is Gaussian), with the 
result 

Z{m,n,c,N) = J f[ dB^ det [I <g> I - c^BW * ® B^]~ n/2 e~ NTr ^ , 

a=l a 

(5.13) 



6 The reader must distinguish between these loops, made of double-lines of a definite color, 
from the oriented loops along which the matrix indices run. 
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where I stands for the N x N identity matrix, Cg> denotes the usual tensor product of 
matrices, and the superscript t stands for the usual matrix transposition. The prefactor 
On is fixed by the condition Z(m, n, c = 0, N) = 1. With this form, it is easy to take the 
logarithm and to let n tend to 0, with the result 



1 f m 

f(m,c,N) = -— / Y[dB (a) Tr(Log[l®I-cJ2B (a) 1 ® B (a) ])e~ NTr ^ 

2 ° N J a =l a 



( B (")) 2 
2 



oo 

C 



P =i 



a=l 

2, 



£ (\^(b^...b^)\ 2 ) 

p=l l<Qi,...,a p <m 



Gauss 



(5.14) 

where we still use the notation (...) Gauss for the multi-Gaussian average over the matrices 
B^ a \ a = l,2,...,m. The modulus square simply comes from the hermiticity of the 
matrices B^ a \ namely 

Tr(JJs( a ^*) = Tr(JjB^*) = Tr(JjB^))* . (5.15) 

Taking the large iV limit ( |5.12|) , it is a known fact |5|] that correlations should factorize, 
namely 

(\Tr(f[B^)\ 2 ) Gauss — |Wns(-))) Gauss | 2 . (5.16) 

1 = 1 1 = 1 

By parity, we see that only even p's give non-vanishing contributions, and comparing with 
( |5.12|) we find a closed expression for the meander numbers of order n with k connected 
components 

n 1 2n 

$>/«m fc = Yl l^^^n^^GaussI 2 . (5.17) 

fc = l l<ai,...,Ct2n<"i i=l 

This expression is only valid for integer values of m, but as it is a polynomial of degree n in 
m (with vanishing constant coefficient), the n first values m = l,2,...,nofm determine it 
completely. So we only have to evaluate the rhs of ( |5.17| ) for these values of m to determine 
all the coefficients Mn ■ 




Fig. 16: The connected toric meander of order 1: it has only 1 bridge. 
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The relation ( |5.17| ) suggests to introduce higher genus meander numbers, denoted 
by M^ k \h), with M 2 ( ^(0) = (note that the indexation is now by the number of 

intersections, or bridges), through the generating function 

oo oo p 

J2J2 M i k) ^ mkN2 ~ 2h = E (|Tr(n5 (ai) )| 2 >Gauss , (5.18) 

h=0 k=l l<oi,...,a p <m i=l 

which incorporates the contribution of all genera in the Gaussian averages. Note that the 
genus h is that of the corresponding black and white graph and not that of the river or 
the road alone. In particular, the river (resp. the road) may be contractible or not in 
meanders of genus h > 0. As an example the = 1 toric meander is represented in 

Fig.16. 




Fig. 17: A typical graph in the computation of the rhs of ( |5.18|) . The two 
p-valent vertices corresponding to the two traces of words are represented 
as racks of p double legs (p = 10 here). The connected components of the 
resulting meander (of genus h = on the example displayed here) correspond 
to loops of matrices B^ a \ This is indicated by a different coloring of the 
various connected components. Summing over all values of ctj yields a factor 
m per connected component, hence m 3 here. 

The relation ( |5.18|) can also be proved directly as follows. Its rhs is a sum over 
correlation functions of the traces of certain words (products of matrices) with themselves. 
More precisely, using the hermiticity of the matrices B^ a \ the complex conjugate of the 
trace Tr(]7 1<i<2n B^ ai ^) can be rewritten as 

Tr( Y[ S (Ql) )* = Tr( JJ fl(«a»+i-0) ? ( 519 ) 

l<i<2n l<i<2n 

i.e. in the form of an analogous trace, with the order of the S's reversed. According 
to the Feynman rules of the previous section in the case of only black matrices, such a 
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correlation can be computed graphically as follows. The two traces correspond to two 
p-valent vertices, and the Gaussian average is computed by summing over all the graphs 
obtained by connecting pairs of legs (themselves made of pairs of oriented double-lines) by 
means of edges. Re-drawing these vertices as small racks of p legs as in Fig. 17, we get a 
sum over all multi-component, multi-genera meanders (compare Fig. 17 with Fig. 5). More 
precisely, the edges can only connect two legs with the same matrix label a, which can be 
interpreted as a color: indeed, we have to sum over all colorings of the graph by means of 
m colors. But this coloring is constrained by the fact that the colors of the legs of the two 
racks have to be identified two by two (the color of both first legs is «i,...,of both p-th legs 
is otp). This means that each connected component of the resulting meander is painted 
with a color a G {1,2, ...,m}. A graph of genus h comes with the usual weight JSf 2 ~ 2h . 
Summing over all the indices ati, ct v = 1, 2, m, we get an extra factor of m for each 
connected component of the corresponding meander, which proves the relation (|5.18j ). 

In the genus case, we must only consider planar graphs, which correspond to genus 
meanders by the above interpretation. Due to the planarity of the graph, the two racks 
of p — 2n legs each are connected to themselves through n edges each, and are no longer 
connected to each other: they form two disjoint arch configurations of order n. This 
explains the factorization mentioned in eq. (|5.16|) , and shows that the genus meanders 
are obtained by the superimposition of two arch configurations. The beauty of eq.( |5.17[ ) 



is precisely to keep track of the number of connected components k in this picture, by the 
m-coloring of the connected components. 

This last interpretation leads to a straightforward generalization of (|5.18|) to semi- 
meanders and many meander-related numbers. 

5.4- Matrix expressions for semi-meanders and more 



In view of the above interpretation, we immediately get the generalization of eq. (|5.17|) 
to semi-meanders as 

n 

k=l l<ax,...,a n <m (5.20) 

N— >oo N 

To get this expression, we have used the m-coloring of the matrices to produce the correct 
rainbow-type connections between the loops of matrices. 
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More generally, this provides a number of matrix integral identities for the following 
generalized semi-meanders. Let us label the arch configurations of order n by a permu- 
tation fx G S2n, the symmetric group over 2n objects, in such a way that if we label 
the bridges of the arch configuration 1,2, ...,2n, the permutation \i indicates the pairs of 
bridges linked by arches, namely, for any i = 1, 2, ...,2n, fx(i) is the bridge linked to % by 
an arch. By definition, \x is made of n cycles of length 2, it is therefore an element of 
the class [2 n ] of S^n- Note that an element of this class generally does not lead to an 
arch configuration, because the most general pairing of bridges has intersecting arches. A 
permutation \i G [2 n ] will be called admissible if it leads to an arch configuration. Let 
Afj, be the arch configuration associated to some admissible \i G [2 n ]. We can define some 
generalized semi-meander number Mn (Afj) associated to A^ as the number of meanders 
of order n with k connected components whose lower arch configuration is A^. With this 
definition, 

Mi k \K n ) = , (5.21) 

where lZ n is the rainbow configuration of order n, associated to the permutation = 
2n + 1 — i, for i = 1, 2, 2n. In other words, M n (A^) is the number of closures by some 
arch configurations of order n of the lower arch configuration A^ which have k connected 
components. Eq. (|5.20| ) extends immediately to 



k = l l<a 1 ,...,a 2n <™ /_ or) N 

a i= a l j.(i)' « = !.■■■. 2" {O.ZZ) 

lim ^-(Tt(B^B^\..B^)) Gauss , 

N— >oo iv 

where the structure of the lower arch configuration A^ is encoded in the conditions ct^ = 
cv^(i) , i = l, 2n, which identifies the colors of the arches according to A^. Some of these 
numbers will be computed in sect. 6. 3 below. 

Higher genus generalizations are straightforward, by simply removing the large N 
limit in the above expressions, namely 



£ MS(A„h)m k N^ h = £ 



fc>l,/l>0 l<"l,-,"2n< m /r no\ 

( Tr(jB (a 1 ) jB (a 2 )_ B (a 2 „) ))Gauss 

Note that the number 2n of bridges is even here, as we are considering arbitrary genus 
closures of a given arch configuration of order n. In genus h = 0, we recover the numbers 
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M 2 ( n } (A^°) = Mi k) (A^) defined by ( Hp . In the particular case A^ = TZ n , this defines 
higher genus semi-meander numbers M^(h) = (1Z n , h), with the correspondence 
M 2 ( n (°) = Mi h) . The resulting higher genus semi-meanders are obtained generically by 
allowing the 4 ends of two given arches to alternate along the river. As the arches cannot 
intersect each other, this requires increasing the genus of the graph. On the other hand, the 
lower rainbow configuration is contractible, hence the genus is also that of the contracted 
graph obtained by the folding process of Fig. 8. 

It is instructive to calculate the sum over all genera of these numbers, while keeping 
track of the numbers of connected components. We simply take N = 1, in which case the 
Gaussian average becomes an ordinary Gaussian average over real scalars (b^\ ... 7 b^ m ') G 
IR m 



m 

M^(A,Mm k = <C>> W )T>Ga U ss 

k>l,h>0 i=l 

30 „2 

,m— 1 „2n 



Jo 

= m(m + 2)(m + 4)...(m + 2n - 2) , 

where the normalization constant has been fixed by the n = 1 case (the result is m). For 
m = 1, the above simply counts the total number of pairings between 2n legs, namely 
(2n — 1)!! = 1.3.5. ..(2n — 1). Note that the result ( |5.24j ) is independent of the lower arch 
configuration, it holds in particular for semi-meanders. The result ( |5.24j ) is a polynomial 
of degree n, with leading coefficient 1 corresponding to the only (genus 0) meander with 
n connected components, obtained by reflecting the lower arch configuration A^ wrt the 
river. For meanders, we simply get 

M^(h)m k = <($>«) 2 ) P ) Gauss 
k>i,h>o i=i (5.25) 

= m(m + 2)(m + A)...(m + 2p-2) . 

Note that in this case the polynomial is of degree p, with leading coefficient 1, corresponding 
to the (genus 1) meander made of a collection of p loops intersecting the river only once. 

5. 5. Computing averages of traces of words in matrix models 



As a warming up, let us first compute the rhs of eq. (|5.20| ) in the case of one matrix 
m = 1, namely 

7n = lim -|-(Tr(£r)> Gauss . (5.26) 
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By parity, we see that 72 S +i = for all integer s. A simple method usually applied for 
computing Gaussian averages uses the so-called loop equations of the matrix model. In 
the case of one matrix, these are obtained as follows. We write that the matrix integral of 
a total derivative vanishes, namely 



- / (IB 



dB 3l L 



= (-N(B s+1 ) kl B tJ + Y,(B r ) kj {B s - r ) u ) G , 



(5.27) 



r=0 

Taking i = I and j = k, and summing over i, j = 1, N, we finally get 



( Tr(B s + 2 ) ) Gauss = 1 J2( Tr(B r ) Tr(B s ~ r ) ) Gauss . (5.28) 

In the large iV limit, due to the abovementioned factorization property, only the even 
powers of B contribute by parity, and setting s + 2 = 2n, this becomes 

n-l 

72n = ^ 72r72n-2r-2 , (5.29) 
r=0 



valid for n > 1, and 70 = 1. This is exactly the defining recursion (|3.1| ) for the Catalan 
numbers, hence 72 n = c n , whereas 72n+i = 0. In this case, the equation ( |5.20|) reduces 
therefore to the sum rule ( |3.8j ) for semi-meanders. Similarly, when m = 1, the equation 
( |5.17|) reduces to the sum rule ( p.7|) for meanders. An important remark is in order. It 
sould not be surprising that the general recursion principle for arches of Fig. 6 resembles 
the large limit of loop equations in matrix models. The role of the leftmost arch is 
played in the latter case by the differentiation wrt the matrix element Bji\ it can act at all 
the matrix positions in the word, which cut it into two even words, and can be graphically 
interpreted as just the right bridge position of the leftmost exterior arch in the pairings of 
matrices necessary to compute the Gaussian average of the trace of the word. 

More generally, the loop equations for the Gaussian m-matrix model enable us to 
derive a general recursion relation for traces of words. The most general average of trace 
of word in m matrices in the large N limit is denoted by 

lpi,P2,---,Pmk 

1 (5.30) 
lim — (Tr(( J B( 1 )) Pl ( J B( 2 )) P2 ...(S( m )) Pm (5 (1) ) Pm+1 ---(5 (m) ) Pmfe ))Gauss • 
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In the above, some powers pj may be zero, but no m consecutive of them vanish (otherwise 
the word could be reduced by erasing the m corresponding pieces). Of course 2p = Y^iPi 
has to be an even number for ( |5.30| ) to be non-zero, by the usual parity argument. For 
m = 1, we recover 7^ = 7 P . If u = exp(2z7r/m) denotes the primitive m-th root of unity, 
then we have the following recursion relation between large iV averages of traces of words, 
for m > 2 

mk — l 

lpi,P2,---,Pmh / 4 lPl,---,Pj 'Pj + l,---,Pm.k ' V / 

When j is not a multiple of m, it is understood in the above that the multiplets (pi, ■■■,Pj) 
and (pj+i, ■■■,Pkm) have to be completed by zeros so a to form sequences of m-uplets. 
For instance, we write 7g 3 ' ) = = 7q 3 q 3 . Note also that if only q < m matrices 
are actually used to write a word, the corresponding 7^ can be reduced to a by 
erasing the spurious zeros (for instance, 7g 3 Q = 'J^)- Together with the initial condition 

7o™..,o,2rH-i,o,---,o = 72n+i = and 7o^ ) ,o,2n,o,...,o = 72n = c n , this gives a compact recursive 
algorithm to compute all the large iV averages of traces of words in any multi-Gaussian 
matrix model. 

A direct proof of eq.( |5.31|) goes as follows. Throughout this argument, we refer to the 
labels a = 1,2, ...,m of the matrices as colors. The quantity lp?;..,p mh can be expressed 
as the sum over all possible planar pairings of matrices of the same color in the corre- 
sponding word, i.e. the sum over all arch configurations of order p (encoded in admissible 
permutations \i of 1, 2, 2p), preserving the color of the matrices (the matrices sitting at 
positions i and have the same color a). Such a color-preserving arch configuration 



appears exactly once in the lhs of (|5.31|) . We will show the relation (|5.31|) by proving that 
each such term also comes with a coefficient 1 in the rhs of (|5.31|) . Let us evaluate the rhs 



of (|5.31| ) in this language. 

The index j may be viewed as the position of a "separator" , which cuts the color-preserving 
arch configurations into two disconnected pieces. The separator positions are labeled by 
the index j = 1,2,..., mk — 1. For a given color-preserving arch configuration in the lhs of 
( |5.31|) , the only terms of the rhs of ( 5.31 ) contributing to it are those where the separator 
index j takes its values at positions inbetween the exterior arches linking various blocks of 
the same color. These positions will be referred to as available positions. 

Available positions satisfy the two following properties, illustrated in Fig. 18. 

(i) any two successive available positions are labelled by successive integers modulo 
m: if an available separator position sits between two blocks of colors a and a + 1, hence 



42 



A.B C A 



B A C A A A B A C 



1 



5 6 7 8 



Fig. 18: A sample color-preserving pairing of matrices for m = 3 matrices. 
Each block of matrices is denoted by a single letter A, £?, C, according to its 
color 1, 2, 3. The precise pairing of matrices within blocks is not indicated 
for simplicity. The available separator positions are indicated by arrows. One 
checks that (i) the positions are consecutive modulo 3, and (ii) the number 
of available positions is 5 = 2x3 — 1. 

at a position of the form j = sm + a, then the block of color a + 1 is linked to other blocks 
of the same color, and the next available separator position sits between a block of color 
a + 1 and a block of color a + 2, hence at a position of the form j = tm + a + 1. These 
two positions are consecutive modulo m. 

(ii) the first available separator position is j = ml + 1: it sits to the right of the 
first set of related blocks of color 1. The last available separator position is of the form 
j = qm + (m — 1), as it sits to the left of the rightmost set of related blocks of color m. 
Hence, thanks to property (i), the total number of available separator positions is of the 
form rm — 1. 

Consequently, the total contribution of a given color-preserving admissible pairing in 
the rhs of ( |5.31| ) is of the form 



So we have proved that the total contribution of each admissible pairing in the rhs of ( |5.31| ) 
is 1. This completes the proof of (|5.31|) . 



average of the trace of a particular word in m = 3 matrices. We wish to compute 



m 



m—l 




(5.32) 



Let us show explicitly how to use the recursion ( |5.31| ) to compute the thermodynamic 



72,1,2,0,1,0 



hm - (Tr (A 2 BC 2 B)) 



auss 5 



(5.33) 
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where we denote by A, B, C the matrices with respective colors 1, 2, 3. Applying ( |5.31] ) 
we get 



(3) - (1) (2) -2 (2) (2) 

72,1,2,0,1,0 = -.772 72,2 - J 72,l72,l 



(3) (1) • (3) (1) . 2 (3) 
72,1,271 - J72,l,27l -J 72,1,2,0,1,0 



(5.34) 



where j = exp(2z7r/3), and the various numbers m of matrices have been reduced to their 
minimal value. R 
integer s, we get 



minimal value. Regrouping the ^y^i 201 o' s ) an< ^ using the fact that 72^1 = 72s+i = for 



( 1 +i 2 )72 3 i, 2 ,o,i,o = -r- "2.2 

Using the recursion (|5.31|) for m = 2, we get 



(2) (1) (1) / n2 

72,2 = 72 72 = (Cl) • 



(5.35) 



(5.36) 



Finally, we get 



(3) 

72,1,2,0,1,0 



( Cl ) 6 = 1 . (5.37) 
This corresponds to the only way of pairing matrices of the same color in the sequence 




A A B C C B 



5.6. Equivalent coloring problems 







(a) 



(b) 



(O 



(d) 



(e) 



Fig. 19: The 5 semi-meanders of order 3. 



n 



The language of sect. 5. 5 above leads naturally to the introduction of the number ft 
of /c-colored multi-component semi-meanders of order n. By /c-colored, we mean that 
exactly k distinct colors are used, but we do not distinguish between colorings differing 
only by a permutation of the colors. For illustration, let us take n = 3. The 5 semi- 
meanders of order 3 are displayed in Fig. 19, and denoted (a), (6), (c), (<i), (e). Only the 

( 3) 

last one (e) can be colored with k = 3 distinct colors, hence p,^ = 1. With k = 2 colors, 
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(2) 

there are 3 ways of coloring (e), and 1 way of coloring (c) and (d), hence a total of fi 3 = 5. 
Finally, each connected semi-meander (a) and (b) can be colored in a unique way with 
k = 1 color, hence = 5. 



In terms of these numbers, the rhs of eq. (|5.20| ) is expressed as a sum over the number 



k of colors used, with a total choice of m(m — l)(m — 2)...(m — k + 1) colors among m, 

- (k) 

weighed by the numbers ix n ■ This implies a relation between coloring numbers \i and 
meander numbers M 

n n 

rh n (m) = J2^i k) m k = ^ ffi m(m - l)...(m - k + 1) . (5.38) 

k=l k=l 

Remarkably, the numbers ft do not carry direct information about the numbers of connected 
components of the semi-meanders, but their knowledge is sufficient to know the M exactly, 
by use of the Stirling numbers. Note also that eq.( |5.38|) expresses the polynomial m n (x) 
in two different bases, {x k , k = 1, n} and {x(x — l)...(x — k + 1), k = 1, n}, of the 
space of degree n polynomials with vanishing constant coefficient. For n = 3, we simply 
have 

x 3 + 2x 2 + 2x = x(x - l)(x - 2) + 5x(x - 1) + 5x . (5.39) 
The same type of expression holds for meanders, namely 

n n 

J2 M n ] ™ k = m(m-l)...(m-k + l) , (5.40) 

k=l k=l 

(k) 

where jXn denotes the number of /c-colored multi-component meanders of order n. 
5.7. Combinatorial expression using symmetric groups 




Fig. 20: An arch configuration of order 3 and the corresponding interpreta- 
tion as a ribbon graph, with V = 1 six-valent vertex and E = 3 edges. On 
the intermediate diagram, the arches have been doubled and oriented. These 
oriented arches indicate the pairing of bridges, i.e. represent the action of 
[i. Similarly, the oriented horizontal segments indicate the action of the shift 
permutation a. Each oriented loop corresponds to a cycle of the permutation 

0[1. 
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Admissibility condition. In sect. 5. 4 above, we have seen how an arch configuration 
of order n could be encoded in an admissible permutation \i belonging to the class [2 n ] 
of Sm- Let us write the admissibility condition explicitly. This condition states that 
arches do not intersect each other, namely that the ribbon graph (see Fig. 20) with only 
one 2n-valent vertex (the 2n bridges), whose legs are connected according to the arch 
configuration, is planar, i.e. of genus h = 0. This graph has V = 1 vertex, and E = n 
edges (arches) . Let us compute its number L of oriented loops in terms of the permutation 
fj,. Let a denote the "shift" permutation, namely a(i) = i + 1, i = 1, 2, 2n — 1 and 
a(2n) = 1. Then an oriented loop in the ribbon graph is readily seen to correspond to 
a cycle of the permutation a[i. Indeed, the total number of loops is L = cycles(cx//), the 
number of cycles of the permutation o[i. The admissibility condition reads 

x = 2 = L - E + V = 1 - n + cycles(a^) 

(5.41) 

cycles(a / u) = n + 1 . 

Note that if we demand that the ribbon graph be of genus h, the above condition becomes 

cycles(a^) = n + 1 - 2h . (5.42) 

Connected components. Given an admissible permutation \i G [2 n ], let us now 
count the number of connected components of the corresponding semi-meander of order 
n. Let r be the "rainbow" permutation r(z) = 2n + 1 — i. Note that r changes the parity of 
the bridge label. On the other hand, the admissible permutation \i is readily seen to also 
change the parity of the bridge labels. As a consequence, the permutation t\i preserves 
the parity of bridge labels. In other words, even bridges are never mixed with odd ones. 
The successive iterations of the permutation t\x describe its cycles. The corresponding 
meander will be connected iff these cycles are maximal, namely fi has two cycles of length 
n (one for even bridges, one for odd bridges), i.e. T/U G [n 2 ]. We get a purely combinatorial 
expression for connected semi-meander numbers 

M n = card{^i G [2 n ] \ cycles(a^) =n+l, and r[i G [n 2 ]} . (5.43) 

More generally, the semi-meander corresponding to [i will have k connected compo- 
nents iff r\i has exactly k pairs of cycles of equal length (one over even bridges, one over 
odd ones). 

46 



Character expressions. The above conditions on various permutations are best 
expressed in terms of the characters of the symmetric group. Denoting by [i Ui ] the class of 
permutations with Vi cycles of length i, and labelling the representations of S^n by Young 
tableaux Y with 2n boxes as customary, the characters can be expressed as 



XY (\i Ui }) = det( Pi+£iH (d.)) 



(5.44) 



where the Young tableau has ii boxes in its z-th line, counted from the top, t v = ]T„ -fr, 
p m (9.) is the m-th Schur polynomial of the variables 6>i, 62,- ■■ 

pm(e.) = £ nfi> ( 5 - 45 ) 

fc i >0,i = l,2,... i *" 

and we used the symbol f{0)\t v for the coefficient of the monomial Y\i -^-\ lIV the polynomial 
/(#.). As group characters, the xy's satisfy the orthogonality relation 

Z>([A])xy(M) = J^SWM , (5-46) 

where the sum extends over all Young tableaux with 2n boxes, [A] denotes the class of a 
permutation A G £2™, and |[A]| the order of the class. The order of the class [i Vi ] is simply 

1™ = ■ (5 - 47) 

The orthogonality relation ( |5.46| ) provides us with a means of expressing any condition on 
classes of permutations in terms of characters. It leads to the following compact expression 
for the connected semi-meander numbers 



r/i] , [n 2 

SA i =n + l 

»-Aii 1 ir^2 



v v v I [211 |^]| IK] I (548) 

((2n)!) 3 

es 2 „ /t*eS 2 « Y,Y',Y" U ; ; 



x Xr(M)xy([2 n ])xy (M)xy(^ l ])xy»([^])xY»([n 2 ]) . 

Analogous expressions hold for (higher genus) semi-meanders with k connected com- 
ponents and for meanders as well. These make completely explicit the calculation of the 
various meander-related numbers. Unfortunately, the characters of the symmetric group 
are not so easy to deal with, and we were not able to use them in an efficient way to 
enhance our numerical data. 
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6. Properties 

In this section, we derive a number of identities and inequalities involving various 
meander and semi-meander numbers. We also compute exactly the number of closings of 
some particular lower arch configurations, made of repetitions of simple motives. 

6.1. Sum rules 



In addition to the sum rules ( p.7|) (|3.8|) and (|5.24|) (|5.25| ) for (higher genus) meander 



and semi-meander numbers, let us derive the following 

n \ j y if n = 2p + 1 

n 



k=i 



(6.1) 

JO if n is even 
c p if n = 2p + 1 



The reader can easily check these sum rules with the data of tables I and II above. 

To prove both sum rules, we introduce the notion of signature of an arch configu- 
ration. Let us denote by |.A| the order of an arch configuration A. The signature of A is 
defined recursively by the initial value on the vacuous configuration and its behaviour under 
the processes (I) and (II) of sect. 4.1, which enable us to construct any arch configuration. 
We set 

sig(0) = 1 

sig((/M) = (-l)W sig(.A) (6.2) 

sig((II).A) = (-1)^ +1 sig(^t) . 

It is important to notice that the signature behaves differently under (I) and (II), hence 
depends on the number of connected components in semi-meanders for instance. More 
precisely, if M. is a semi-meander of order n with k connected components, Ai e SM n ^, 
and V € SM n ^ 1 then they both originate from the single arch of order 1 (with signature 
— 1) by a total of (n — 1) actions of (I) and (II), each adding one arch, except that V is 
obtained by acting (n — 1) times with (I), whereas M. is obtained by acting (k — 1) times 
with (II), and (n — k) times with (I). Due to ( |6.2|) , this implies a difference in the signatures 

sig(M) = (-1)*" 1 sig(P) , (6.3) 

whereas after (n — 1) actions of (I), adding one arch each time, the signature of V reads 

-, , n ( n — 1 ) 

sig(P) = (-1) 1+ — ^ . (6.4) 
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Thus, the signature of a semi-meander only depends on its order n and its number of 
connected components k, namely 

sig(M) = (-1)"+^, VMeSM n , k . (6.5) 

In fact this signature can also be easily related to the numbers of arches of given depth 
in A as follows 

sig(^) = (_l)£;>i iA ^ , (6.6) 
where A(A, j) (see sect. 3) denotes the number of arches of depth j in A. One readily 



checks that eq.( B.6|) solves the recursion relations (|6.2|) . For instance, the signature of the 
rainbow configuration of order n is 

sig(^) = (-1)^ , (6.7) 

in agreement with (|6.3|) for k = n. To conclude, the signature of a semi-meander only 
depends on its order n and its number of connected components k, namely 

sig(TW) = (-l)^ 2 ^, \/MeSM n , k . (6.8) 

The sum rule ( |6.1| ) for semi-meanders is obtained by computing the sum of all the 
signatures of multi-component semi-meanders of order n. Let us introduce 



s(n) = E sig (^) 
AeA n 

fc=i 



AeA n 

(6.9) 



where the first sum extends over all the arch configurations of order n. Using the usual 
recursive reasoning of Fig. 6, the leftmost arch of a configuration separates the configuration 
into two configurations of smaller orders, say X of order j and y of order n — 1 — j (the 
only configuration of order is just 0) . As the j arches of X are shifted downwards in the 
global configuration, we obtain the recursion relation 

n-1 

s(n) = E E E sig(^)sig(y)(-l) J+1 

j=o xeA 3 yeA B _i_j 

= J2(-iy+h(j) s (n-l-j). 

j=o 

49 



With s(0) = 1, we immediately find that s{2p) = for p = 1,2, and that 

p-i 

s(2p + l) = £s(2j + l) s(2p-2j-l) , (6.11) 

j=0 

with s(l) = —1, hence for n = 2p + 1 

s(2p + l) = (-l) p+1 c p . (6.12) 

Comparing with (|6.9|), and noticing that 1 + n(n — l)/2 = 1 + (2p + l)(2p)/2 = p + 1 
modulo 2, this completes the proof of ( |6.1| ) for semi-meanders. 

In the case of meanders, let us define the signature of a given meander Ai with upper 
and lower arch configurations *A up and ^4down respectively, as the product^ 

sig(M) = sig(Ap) sig(^l down ) . (6.13) 

With this definition, we find that the signature of a meander only depends on its order 
and on its number of connected components through 

s i g (M) = (-l) k+n . (6.14) 



Therefore, 

n 

£ M W (_!)*"! = (_!)-! E ^ig(Ap) Sig(^down) 

= (-l)"" 1 ^) 2 ( 6 - 15 ) 

f if n = 2p 

\(c p ) 2 ifn = 2p+l 

This completes the proof of ( |6.1|) for meanders. 

7 Note that in a previous definition, the signature of a semi-meander was referring to its upper 
arch configuration only. Considered as a particular meander, its signature is simply multiplied by 
that of the raibow of order n ( |6.7| ) . 
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Combining (|3.7|) and (|3.8j ) with (|TT])> we get the equivalent set of rules 



n 2 n 

EC = f = 5>£" 

k=l k=l 

n+1 2 i 2 

E M (2k-l) _ C 2n +1 + C n 

IV1 2n+l — 2 

k=l 

n 2 2 

E A# (2fc) _ C 2n+1 - C n 

IV1 2n+l — 2 



(6.16) 



fc=i fc=i 

n+l 

E,->(2A;-1) _ C2n+1 T C n 
iW 2n+l - 2^ 

fe=l 

(2fc) _ C2n+1 — C n 



fc=i 

tf.i?. Inequalities 

The above sum rules ( |6.16| ) lead to some immediate inequalities, by simply using the 
positivity of all the meandric numbers involved, namely 

2n - 2 

M (2k-1) < C 2n+1 + g 
Zn-\- 1 — 2 

2 2 

2 - +1 S 2 (6.17) 

(fc) . c 2n 



-^2n < 



2 



,-v(2fe-l) ^ C2n+1 + C n 
M 2n+1 — 2 

»>(2fc) ^ C 2n +1 - C n 
iW 2n+l — 2 ' 

for all allowed fc's. 

In particular, according to the leading behaviour c n ~ 4 n when n — > oo, this gives an 
upper bound on the numbers R and .R ( |2.12|) (|2.13|) 



R < 16 i? < 4 . (6.18) 
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Defining as before Rk and Rk as the leading terms in respectively meander and semi- 
meander numbers with k connected components for large order n (with R\ = R and 
Ri = R) 

„ (i2 fc )» M( fc ) ~ (R k y , (6.19) 

this also gives the bounds 

Rk < 16 R k < 4 . (6.20) 

The matrix model formulation of sect. 5 gives rise to more interesting inequalities, 
relating meander and semi-meander numbers. They are based on the Cauchy-Schwarz 
inequality for the hermitian product (A, B) over complex matrices defined by 

(A,B) = Tr(AB f ) , (6.21) 

where the superscript f stands for the hermitian conjugate = B 1 * , namely 

\(A B)\ 2 < (A, A) (B,B) . (6.22) 

Starting from the expression ( |5.17|) , let us write 

A = S (ai) B (a 2 )__ B (ai) &nd g\ = fl (a,+i) .. -jB (a2„) ? (g_ 2 3) 

then ( |6.22| ) becomes 

2n 11 2n 2n-l 

|Tr(JJS (Ql) )| 2 < Tr( JJs( ai) JJs( ai + 1 -^ ) )Tr( j f B^ \\ fl(«a»+i-i)J . ( 6 . 2 4) 
i=l i=l j=l i=/+l i=l 

Taking the thermodynamic Gaussian averages of both sides and summing over the colors 
cci = 1, 2, m, we finally get 

n I In— I 

£ M<*> m fc < ( £ m«) ( £ ^ m') , (6.25) 

fc=l i=l j=l 



where the rhs has been identified by use of (|5.20|) . This inequality between polynomials 



results in many inequalities for the meander and semi-meander numbers themselves. 
In the trivial case I — 0, i.e. for the choice A = I, we have 



2n 

E^' "»* < E^2n m k ■ (6-26) 



fc=l k=l 
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This inequality is in fact obviously satisfied coefficient by coefficient for k = 1,2, n, as a 
meander of order n is a particular case of semi-meander of order 2n (which does not wind 
around the source of the river), namely 



M^ k) < M 2 ( " } Vfc = l,2,...,n. 



(6.27) 



This implies the inequalities between the leading terms Q6.2U ) 



R k <R 2 k Vfc = l,2,... 



(6.28) 



Before taking the large N limit, (|6.24| ) also leads to an inequality involving higher 
genus meander and semi-meander numbers 



M^ k \h) m k N 2 ~ 2h < ( m k N 1 ~ 2h ) 

h>0,k>l 



h>0,k>l 



(6.29) 



h>Q,k>l 




Fig. 21: The algorithm for closing a given upper arch configuration into a 
connected meander. The upper arch configuration A is of order 6. The path 
starts on the leftmost exterior arch of A: the first lower arch connects it to 
the second exterior arch, and is numbered f . Then the path goes to the left, 
and connects the depth 2 arches of A through lower arches number 2, 3 and 
4. The leftmost depth 2 arch of A is connected to the depth 3 one through 
the lower arch number 5, and finally the depth 3 arch of A is connected to 
the leftmost exterior arch of A through the lower arch number 6. 
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A minoration of the connected meander numbers M n is obtained as follows. For any 
given upper arch configuration A, we have a simple algorithm for constructing its closing 
by a lower arch configuration, such that the resulting meander is connected. Consider the 
exterior arches of A, and draw lower arches connecting any two consecutive exterior arches, 
as shown in Fig. 21: this creates a path (oriented, say from the left to the right) connecting 
the exterior arches of A. In a second step, the path is continued by connecting from the 
right to the left the arches of depth 2 in A, through lower arches which are either disjoint 
from the previous ones, or contain some of them. The lower arch configuration is further 
continued by connecting the arches of increasing depth j in A, through lower arches. The 
resulting path zig-zags through the arches of A, describing those of even depth from right 
to left and those of odd depth from left to right. Once the last arch of A is described, a last 
lower arch is constructed, connecting this last arch to the leftmost exterior arch we started 
with. The total being closed and going through all the arches of A as well as the lower 
ones, we have constructed a connected meander. Therefore the total number of connected 
meanders of order n, M n , is larger than that of arch configurations of order n, c n 

M n > c n . (6.30) 

This implies that 

R > 4 R > 2 , (6.31) 
where we have used the inequality (|6.28|) for k = 1. 

6.3. Special cases 

In this section, we compute exactly some of the generalized semi-meander numbers 
of sect. 5. 4. We will consider the closings of some particular lower arch configurations 
Afj,. It turns out that the numbers can always be found by relatively simple recursions 
provided the arch configuration A^ is the repetition (say p times) of some small order arch 
configuration. 

Single arch. The simplest example is Ai, obtained as the repetition (n times) of single 
arches of order 1, and corresponding to the permutation \i: /x(2i — 1) = 2i, % = 1, 2, n. 
Let us compute the corresponding numbers M^\Ai) by the method of loop equations 
described in sect. 5. 5. The generating function a n (m) for these numbers satisfies ( |5.22| ), 
namely 

n 1 m 

a n (m) = J2Mi k) (A l )m k = jim -<Tr( ^(B«)T>g mss • (6.32) 

k=l j = l 
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Let us use the loop equation method of sect. 5. 4 to derive a recursion relation for a n (m), 
by writing that the multi-matrix integral of a total derivative vanishes 



/TO rj m . . 

II dB ' a) 7^m [[£( B(r, ) 2 r lBW U- NTrE "' * . (6-33) 

a=l ® ji r=l 

hence, taking i = I and j = k, and summing over i,j = 1,2, ...,N 

m m 

iV(TV([^( J BW) 2 ]- 1 (S^) 2 )) Gauss = (MHE(B (r) ) 2 ] B_1 ) 

n — 1 

+ Tr([^(SM) 2 p- 1 S^)Tr([^(SW) 2 ] n - 1 ^S^) (6 " 34) 

+ TV([^( J BM) 2 F- 1 )TV([^( J BW) 2 ]- 1 ^( J B^) 2 )) Gauss . 

By the factorization property ( |5.16| ), the term on the second line vanishes by imparity. 
Summing then over (3 = 1, 2, m and letting iV —>■ oo, we finally get a recursion relation 
for a n (m) ( |6.32| ) 

n-l 

ot n {m) = (m - l)o! w _i(m) + aj(m) a n _!_j(m) , (6.35) 

i=o 

whith the initial condition «o(w) = 1. This determines uniquely the polynomials a n (m) 
of degree n in m. Remarkably, we find 

n 

a n (m) = ^2l(n,k) m k , (6.36) 

k=i 

where the numbers I(n, k) are given by ( |3.16| ) ! In particular, we find the average number 
of connected components in closings of Ai 

(compi) „ = T.UkI(n,k) = n + 1 (6 37) 

The result (|6.36| ) is readily checked by multiplying both sides of the recursion relation ( |3.15| ) 
by m k and summing over k, which yields ( |6.35| ). So the number of arch configurations of 
order n with k interior arches is also the number of closings of A\ with k connected 
components. Note that there is J(n, 1) = 1 connected closing of A\ only, which may be 
obtained by the closing algorithm presented at the end of the previous section. 
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This remarkable coincidence can be explained as follows. Let us constuct a one-to-one 
mapping U between arch configurations, which sends arch configurations which close A± 
with k connected components onto arch configurations with k interior arches. The map U 
is defined recursively as follows 




(6.38) 



Fig. 22: The mapping U between arch configurations is indicated by arrows, 
for the cases of order 1, 2, 3. In these cases U appears to be involutive, but 
this is only an accident, as one readily checks for n > 4. 

A few examples of the action of U are given in Fig. 22. The coincidence between the 
number of connected components after closing a given configuration by A± , and the number 
of interior arches in its image by U needs only to be checked on a configuration with 
one exterior arch, because the second relation of ( |6.38| ) can be used to reduce any arch 
configuration into a succession of arch configurations with one exterior arch only, whereas 
the numbers of connected components in the closure by A% and of interior arches both add 
up in a succession of arch configurations. Counting the number k of connected components 
in the lhs of the last line of ( |6.38| ), we find 

k = #c.c.(/^Vy\) = #c.c.(I)+#c.c.(^) (6.39) 
On the other hand, the number of interior arches in the rhs is 

i = #mt.U(X) + #int.C/(^Y\) . (6.40) 
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Therefore, comparing ( |6.39| ) and Q6.4Q ), we get the desired result by recursion on the order 
of the arch configuration: the number of connected components in the closing by A± is 
equal to the number of interior arches after the action of U. 

Double arch. The repetition (say n times) of a double arch (rainbow of order 2) gives 
an arch configuration Ai of order 2n. Let us compute the number (A2) of closings of 
this configuration with k connected components. As pointed out in the previous section, 
the loop equation method is parallel to the general recursion principle of Fig. 6 for arch 
configurations. To further illustrate this similarity, let us compute (A2) in the lan- 
guage of arches, by deriving some recursion relations. When trying to apply the general 
process of Fig. 6, we immediately see that the recursion will not involve closings of A2 only, 
but requires the introduction of some new arch configuration closings. 
Let us introduce the notations 



A 



v. 



(k) 

2n+l 
(k) 

(k) 
2n+3 



M^ +1 (A 2 

m£Ha 2 ) 



(6.41) 



M. 



2n+3 



for n > 0. The configuration A 2 and its slight modifications listed in (|6.41D form a system 
which is closed under the recursion principle of Fig. 6. 




Fig. 23: The three recursion relations involving the numbers ^n+ii A*2n+2 
(k) 

and z4n + 3 respectively. The leftmost arch is indicated in thick solid line. The 
symbol E indicates a sum over the position of the right bridge of the leftmost 
arch. 
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More precisely, we find the three recursion relations, graphically depicted in Fig. 23 

n— 1 k 

A 2n+1 — P2n P2n "> / , / , A 2i+l A 2n-2i-l 

j=0 i=l 

n-1 fc-1 
3=1 i=l 

(As) _ x (fc) (fc-2) (fc) 

f x 2n+2 — A 2n+1 ^2n + ^2n+l 



(6.42) 



n— 1 fc n fc — 1 

A 2j + 1^2n-2j ^2n+l-2j 

(fc) _ (fc-1) x(fe-l) , (fc) 
^2n+3 ~~ A*2n+2 A 2n+1 + ^2n+2 

n k n k 

, w„w \ (*+*-*) i W // (l) // (fc+1_l) 

"^Z^Zi V 2j+l A 2n+l-2j ^2j ^2n+2-2j ) 

j=0 i=l J = l £=1 

valid for n > and fc > 1. In terms of the generating functions 

2n+l 

x(z, X ) = eea^v** 1 = **+... 

n>0 fc=l 
2n+2 

V(z,x) = E ^ fc n ) +2^ 2n+2 =^(^ + a;2 ) + - (6-43) 

n>0 fc=l 
2n+3 

x) = ]T ]T ^+3^^ 2n+3 = * 3 (* + 3x 2 + x 3 ) + ... 

n>0 fe=l 

we find the following system of algebraic equations 

X(z, x)(- - ___L_1) = (i -f ^(^ + n(z, x)) 

Z X 

( fX (z,x)+x)(--^^-) = (l + fi(z,x))(v(z,x)+zx 2 ) + - (6.44) 
z x z 

(z^,x)+zx 2 )(--^^) = (x + fi(z,x))(x+^^-). 

Combining these three equations, we easily eliminate the functions A and v and get a 
polynomial equation of order 4 for fi(z, x) 

z 2 (l +/i)(x + pb){x + x 2 + 2/j) 2 - xfi[x + x 2 + ft) = , (6.45) 

which, together with the initial condition ii(z, x) = z 2 (x + x 2 ) + ..., determines the numbers 
\x^n completely. The functions A and v are simply expressed in terms of [i thanks to ( |6.44| ) . 
For instance, 

\(z,x) = x ^ x \ ^, (6.46) 

v ' ; z(x + x 2 + 2 / u(z,x)) ' v ! 
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in the sense of formal series expansion in powers of z. 

It is interesting to check the validity of the result ( |6.44[ ) by taking x = 1 , in which case 
the numbers 'Yli<k<2n /"4« = c 2n simply count the total number of arch configurations of 
order 2n, namely 

V(z,l) = J2 c ^z 2n = C+(z)-l, (6.47) 

n>l 

where C_|_ is the even part of the Catalan generating function (|3.3p 

c+w = VTTjWIEg (64g) 

The relation ( |6.44| ) reduces for x = 1 to 

4z 2 C+(z) 4 = C+(z) 2 -l , (6.49) 

in agreement with ( 6.48 ). For x = — 1, we find that fj,(z, — 1) = 0, which translates into 
the sum rule 

2n 

fc=l 

Let us now use the equation ( |6.44|) to derive the leading behavior of the average 
number of connected components over closings of A<i- This number reads 

<comp 2 ) 2ri = ^ = ^ ^ . (6.51) 

We extract from (|6.44|) the generating function fi(z) = d x fi(z, x)\ x —i of fi2n by differenti- 
ating it wrt x and then setting x = 1 

/|M (C + (^)+4)(C + (^)-l)-16^C + (z) 3 

From the study of the limit z —>■ l/4~ in this last expression, we find the asymptotics of 
fi2n namely 

H(z) 



2Jl - 4z 2 . 

I K ' 6.53 
V2-1 4 2n 

* ^ ~ V2 v^(2n)i • 

Dividing this by the large n asymptotics of C2 n ~ 4 2n /( v /7r(2n) 3 / 2 ), we finally get the 

leading behavior of the average number of connected components in closings of A2 

(comp 2 ) 2n = 2n(l - -=) . (6.54) 
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This is to be compared with the analogous expression ( |6.37|) for A\ closings, which gives 
asymptotically (comp 1 ) n ~ n/2. There are much less connected components per bridge in 
average closings of A 2 than A\. 

The connected numbers /4n = M 2n (A 2 ) are generated by the coefficient of x in the 
series expansion of /i(z, x), namely 

00 

fi(z,x) = xfi(z)+0(x 2 ) = x^M 2n (A 2 ) z 2n + 0(x 2 ) , (6.55) 

n=l 

hence as a consequence of ( |6.45| ), the series /j,(z) satisfies an algebraic equation of order 2 



z 2 (l + 2n(z)) 2 - n(z) = 0, (6.56) 



which we can rewrite as 



2z 2 (2/j(z) + l) 2 = (2/i(z) + 1) - 1 . (6.57) 

We recognize the second degree equation for the generating function for Catalan numbers 
upon the substitution 2/j 1 (z) + 1 = C(2z 2 ). We deduce that 



^ = M 2n (A 2 ) = 2 n - x c n V n > 1 . (6.58) 
Analogously, we find that 

\(!) _ <yn„ 

(6.59) 



\W _ on 
A 2n+1 — z L rt 

^2n+l — 1 \ c n+l 



The , 1 connected closings of A 2 completed by a single arch are obtained, starting from 
the A 2 connected closings of order (2n), by the process (I) of Fig. 10. Indeed, the resulting 
meander is of order (2n + 1) and has a single lower exterior arch surrounding A 2 . By 
cyclicity along the infinite river, this arch can be sent to the right of A 2 , thus forming one 
of the A2^ +1 meanders. Therefore, X^n+i * s ec l ua l to the total number of exterior arches in 
connected closings of A 2 (each exterior arch gives rise to one action by (I)). The average 
number of exterior arches in connected closings of A 2 reads therefore 

A (i) 

(ext 2 ) 2n = = 2 , (6.60) 
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independent of n. On the other hand, due to the asymptotic behaviour of the Catalan 
numbers (c n ~ 4 n , n — > oo), we find that for large n 

M 2n („4 2 ) ~ (2V2) 2n • (6.61) 

More cases. The triple arch case .A3 (repetition n times of the rainbow TZ 3 of order 3) 
involves the introduction of 9 slightly different arch configurations. We will not give the 
details of the solution here, and simply give the result: the generating function for the 
number of connected closings of As, 

6{z) = J^MsniAs) z 3n , (6.62) 

n>l 

satisfies the polynomial equation of order 4 

2^ 3 (1 + 6(z))(l + 39(z)) 3 - 6(z)(l + 29(z)) 2 = . (6.63) 

The asymptotic behaviour of the numbers Ms n {Az) is found by solving the coupled system 
of equations 

2z 3 (l + #)(l + 3#) 3 -#(l + 2#) 2 = 

d ( 6 - 64 ) 

— (2z 3 (l + #)(l + 3#) 3 -#(l + 2#) 2 ) = 0. 
ou 

The solution (9*,z*) is such that 

M 3n (A 3 ) ~ (z*)- 3n = (2. 842923.. ) 3n . (6.65) 

The above behaviour is to be compared with that of connected closings of single arch 
repetitions ~ l n and double arch repetitions ~ (2y / 2) 2n = (2.828..) 2n . We see that very 
little entropy has been gained when replacing double arches by triple arches. Intuitively, 
these numbers should form an increasing sequence (with increasing order of the elementary 
raibow which is repeated), which may very well converge to the number R, governing the 
large n asymptotics of semi-meander numbers. 

Another case of interest is the repetition B p (n times) of the arch configuration of 
order p obtained as a repetition of (p — 1) single arches of order 1 below a single arch, 
i.e. corresponding to the permutation ^(1) = 2p, n(2i) = 2i + 1, i = 1,2, ...,p — 1. The 
generating function 

00 

(3(z) = Y,M np (B p ) z np (6.66) 

n=l 
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for the numbers of connected closings of B p can be shown to satisfy a polynomial equation 
of order p 

z p (l + p[3(z)) p -P(z) = 0. (6.67) 
This leads to the large n behaviour 

M np (B p ) ~ (p E ^(p-l)^ E ) np . (6.68) 

For p = 2, we recover the result ( |6.61| ) for double arch repetitions. 

Note that, in view of the above cases, an empiric rule seems to emerge, indicating 
that, for fixed order, the number of connected closings of a motive is increased by nesting 
the arches rather than juxtaposing them. If this is true, it means that the number of 
closings of a given motive is maximal with the rainbow (semi-meander) . We have checked 
this conjecture numerically up to order 12. 

7. Irreducible meanders, exact results 

A multi-component meander is said to be k— reducible if one (proper) subset of its 
connected components can be detached from it by cutting the river k times between the 
bridges. A multi-component meander is said to be k— irreducible if it is not /c-reducible, 
i.e. if no (proper) subset of its connected components can be detached from it by cutting 
the river k times between the bridges. 




(a) (b) (c) 



Fig. 24: Reducibility and irreducibility for meanders. The cuts reducing the 
meanders are indicated by arrows. The meander (a) is 1-reducible, i.e. the 
succession of two meanders along the river. The meander (b) is 1-irreducible 
and 2-reducible. The meander (c) is 1-and 2-irreducible and 3-reducible. 

In Fig. 24, we give a few examples to illustrate the notion of reducibility and irreducibility 
of meanders. The same definition applies to the semi-meanders in the formulation with a 
semi-infinite river. 
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7.1. 1 -irreducible meanders and semi-meanders 



A meander is 1-irreducible if it is not the succession along the river of at least two 
meanders. We can enumerate all the meanders by their growing number of 1-irreducible 

(k) 

components. Denoting by PA the total number of 1-irreducible meanders of order n with 
k connected components, we compute 



7i^+n2 + ...~ra 

fci+fea+. ..=fe 



hence the generating functions 



M(Z,X) = ^ 



z 2n x k 



n>k>l 



(k) 2n k 

n A 



M(z,x) = ^ ] 



z 2n x k 



n>k>l 



P(z,x) = P , 



(k) 2n k 
n * ^ 



n>k>l 



(7.2) 



P(z,x) = Yl P i 

n>k>l 

satisfy the equation 

M(z,x) = P(z,x) + P(z,x) 2 + P(z,xf + ... = - P ^\ X) . . (7.3) 

1 — P(z, X) 

An analogous reasoning for semi-meanders leads to the relation 

P(z,x) 
1 - P(z,x) 

between the generating functions 



M(z,x), (7.4) 



(7.5) 



of respectively semi-meander and 1-irreducible semi-meander numbers of order n with k 
connected components, and P(z, x) defined in (|7.2Q . 

7.2. 2-irreducible meanders and semi-meanders 

The 2-irreducible meanders have been studied in [|1| extensively (under the name of 
irreducible meanders). They use the following equivalent characterization of 2-irreducible 
meanders: for any subset of connected components of a 2-irreducible meander, its set of 
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bridges is not consecutive. Indeed, if it were not the case, one could cut the river before the 
first bridge and after the last one and detach the corresponding piece. One easily checks 
on Fig. 24 (c) that the two sets of bridges of the two connected components of the meander 
are intertwined. 

For completeness, we reproduce here their computation of the numbers q n of 2- 
irreducible meanders of order n, by slightly generalizing their argument to include the 

(k) 

numbers Q n of 2-irreducible meanders of order n with k connected components. The 

(k) 

idea is to enumerate the M n meanders by focussing on their leftmost 2-irreducible com- 
ponent, namely the largest 2-irreducible subset of its connected components, containing 
the leftmost bridge. 




Fig. 25: A general meander. The leftmost 2-irreducible component is rep- 
resented in thick solid lines. The positions marked by # can be decorated 
with any meanders to get the most general meander with this leftmost 2- 
irreducible component. 

In Fig. 25, the leftmost 2-irreducible piece of a meander is depicted in thick solid line. 
Suppose this piece is of order p and has / connected components. The most general meander 
having this leftmost 2-irreducible piece can be obtained by decorating any segment of river 
between two consecutive bridges with arbitrary meanders of respective orders m, n2, ... 
n 2 P , P + n>i + ... + n 2p = n (there are 2p positions, indicated by the symbols # in Fig. 25, 
for these possible decorations), and with respective numbers of connected components k±, 
&2, ••• &2 P , with I + k{ = k. We get the relation 

P+»l +■ ■ ■+"2p= Tl 
l + ki + . . .+k 2p = k 

with the convention that = 1. In terms of the generating functions M(z, x) of eq. 
and 

Q(Z,X) = z2Uxk » ( 7 - 7 ) 

n>fc>l 
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A) 



the relation (|7.6| ) reads 

M(z,x) = Q(z(l + M(z,x)),x) 



In the special = 1, if we denote by 

oo 

B(z) = 1+M(z,l) = ^(c 

n=0 
oo 

q(z) = 1 + Q(z,l) = Y,<ln 

then (17.81) reduces to 



) 2 z 2n 



z 2n 



B{z) = q{zB{z)) . 

The radius of convergence of the series B(z) is 1/4, due to the asymptotics of ( 
that of q(z) is 



n>0 



1 X 1 , 

1 ; (lx -C(-)C{ 1 



4 / 2inx 4 4x 

— (1 - VT^)(1 - v 7 ! - l/x) 
inx 

dx 



2mx 

■ 2tt 



-1 + ^2 - x - l/x) 



1 r 9 

-J o ^ ( -i + 2si„-; 



4 - 7T 
7T 

Hence a leading behaviour for large n [Ij 

2n 



9n 



7T 



4 - 7T / 

This gives an upper bound on the leading behaviour of M n 

U < ( — ^— ^ = 13.3923. 

4 - 7T J 
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(7.14) 



Let us extend these considerations to the semi-meander case. If Q n denotes the 
number of 2-irreducible semi-meanders of order n with k connected components, let us 
enumerate the MA meanders of order n with k connected components, focussing on their 
leftmost 2-irreducible piece, say of order p with I connected components. 



Fig. 26: A general semi-meander. The leftmost 2-irreducible component is 
represented in thick solid lines. The positions marked by can be decorated 
with any meanders to get the most general semi-meander with this leftmost 
2-irreducible component. The last position, indicated by @, can be decorated 
with any semi-meander, possibly winding around the source of the river (*). 

The most general semi-meander with given leftmost 2-irreducible component (as indicated 
by thick solid lines in Fig. 26) is obtained by decorating the segments of river between any 
two consecutive bridges (there are (p— 1) such positions) with meanders of respective orders 
ni, nz,... n p -i, with respectively ki, ...,k p -\ connected components, and the segment of 
river between the source and the first bridge by any semi-meander of order q, with r 
connected components, such that p + q + 2J2 n i = n an d I + r + Y^ki = k. This amounts 
to the following relation 



where the sum extends over non-negative values of rii and ki, with the convention that 
Mq '* = 1 = <3o°' ) = Mq " 1 . This translates into a relation between the generating functions 







(7.15) 



P + <2 + 2(ni + ...+Ti p _i) = Ji 
l+r+ki + . . .+k p _ i —k 




n>k>l 



(7.16) 



Z 



n 



X 



k 



n>k>l 
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and M(z,x) of eq.( [7.2|) (note that in ( [7.16|) we have given a weight z per bridge in the 
semi-infinite river framework, while there would be twice that number of bridges in the 
rainbow closing framework) 



M(z,x) = Q(z(l + M(z,x)),x) \ 
M(z,x 



(7.17) 



■(1 + M(z,x)) = Q(z(l + M(z,x)),x) . 



1 + M(z,x) 

In the special 1, with the generating function 



q(z) = 1 + Q(z,l) = J^Qn z n , (7.18) 



n>0 



where q n = Y^k=i Q™ * s ^ ne total number of 2-irreducible multi-component semi- 
meanders of order n, and with C(z) = 1 + M(z, 1), B(z) = 1 + M(z, 1) defined respectively 
in ( jOD ( |7T9D , the relation ( |7T?D reduces to 

l + zB(z)C(z) = q(zB(z)) , (7.19) 



where we used the fact that (C(z) — 1)/C(z) = zC(z) fl3T4|). Reasoning as above, we find 
the convergence radius of the series q(z), z* = _B(l/4)/4 = (4 — n)/n. Hence we have the 
asymptotics 

(7.20) 



4 - 7T 

This implies the upper bound on the leading behaviour of M, 



R < = 3.659... (7.21) 

4 — 7T 



This upper bound is below the mean field estimate (|4.52| ) which therefore cannot be exact 



and needs to be improved. Again, a more refined study of the relation (|7.19|) gives the 
subleading power law 



~ — 3 



1 ( TV 

3~ 
712 



4 -^y (7.22) 
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8. Conclusion 



In this paper we have studied the statistics of arch configurations, meanders and semi- 
meanders. This study emphasizes the role of exterior arches, instrumental in the recursive 
generation of both semi-meanders and arch configurations, and whose average number 
is directly linked to the corresponding entropies. A complete solution of these problems 
however requires the knowledge of the correlation between these exterior arches and those 
of a given depth (see for instance eqs. (|4.11|) ( |4.12| )). Discarding these correlations leads to 
a rough mean field result. 

Using the alternative formulation (|5.17|) -( |5.20| ) in the framework of random matrix 
models, we are left with the computation of particular Gaussian matrix averages of traces 
of words. Equation (|5.31|) provides a recursive way of computing any such average. As a 
by-product, we have derived an exact formula ( |5.43|) expressing the semi-meander numbers 
in terms of the characters of the symmetric group. 

Finally, we have presented the exact solutions for simpler meandric problems such as 
the closing of particular lower arch configurations (made of rereated motives), as well as 
the 2-irreducible meanders and semi-meanders. 

In conclusion, we are still lacking of an efficient treatment of the main recursion 
relation for semi-meanders. Beyond the entropy problem, it would also be interesting to 
determine values of critical exponents in this problem, such as the exponent governing the 
subleading large order behaviour of the meander and semi-meander numbers ( |2.11| ). The 
latter are nothing but the usual a and 7 configuration exponents of polymer chains. From 
this point of view, a fundamental issue is the determination of the universality class of the 
(self-avoiding) chain folding problem. 
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Appendix A. Combinatorial identities for arch statistics 

Many combinatorial identities can be proved directly by identifying the result as a 
certain coefficient, say of x n , in a (Laurent) series expansion of some function f(x) around 
x = 0. In the following, we use the notation 

'<*>!- s / ^ 
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for the coefficient of x n in the small x expansion of the function /. 

Proof of eq. (|3.12|) . We want to prove the identity (|3.12| ), which may be rewritten as 



E 

k=i 



k 



2n-k 



2n — k\ (k 



n 



n — m fn + m — 1\ (n — m 



n-l 

E 

m=0 

21 + 1 (2n + l 



n 



m 



2n + 1 V n - I 



(A.2) 



where we changed the summation variable to m = n — k. The above sum over m can be 
expressed as 



E 

r, m > 
r-\-m = n — I 



n + m — 1 \ r + I ( r + I 
x 



rn 



n 



(A.3) 



where we introduced a second summation over r = n — I — m. Using the infinite series 
expansion 



1 



E 



m + a — 1 



■/?? 



x 



1 - x) a 

x ' m=U " 

we identify the expression ( |A.3| ) as the coefficient of x n ~ l of some series 

1 



(A.4) 



1 l n d 
■x—(l + x- 



[1 — x) n n dx (1 — x) l+1 
I 



n{l - x) n + l + x 



+ 



l + l 



n(l - x ) n + l + 2 



lfn + l + n-l\ l + lfn + l + l + n-l-1 
H 



it \ n — l ) n V n — I — 1 

I (J(n+J + 1) + (J + 1)(n _ 0) _^_ 

2/ + 1 /2n + 1 



(A.5) 



2n + 1 V n - I 
This proves the identity ( |3.12| ). 

Proof of the identity ( |3.17| ). We want to prove the identity ( |3.17| ), which reads 



n 

n ^ 

k=i 



k\ (n 



I \k \k-l 



n 



1 (n 
nil 



n — I 



E 

k=l 

1 fn\f 2n-l 
nil J In- l + l 



k - I \k - 1 



n 



(A.6) 
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The sum over k in the second line of ( |A.6|) is nothing but the coefficient of x l 1 of some 
Laurent series, identified as 



k=i 



n — l\ ( n 



k- I \k- 1 



(l + x) n "'(l + 



1 



X 



= (l + x) 2n " z | 
This completes the proof of the identity (|3.17|) . 



p 71 — I 



2n - I 
n-l + 1 



(A.7) 
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